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Introduction 


This  paper  documents  a  series  of  subroutines  developed  in 
support  of  inventory  policy  research  concerning  the  impacts  of 
lead  time  variability  upon  D062  Economic  Order  Quantity  Items. 

The  subroutines  are  designed  to  operate  within  the  framework 
of  the  HEDGSIM  Long  Supply  Simulation  Model.  Appendix  A  contains 
Job  Control  Language  (JCL)  statements  required  to  utilize  these 
new  routines  with  other  members  of  the  HEDGSIM  simulation  program 
library.  On  the  other  hand.  Appendix  B  presents  the  Fortran 
source  program  listings  of  new  routines  developed  for  this  study. 

The  new  routines  documented  in  this  paper  may  be  assigned 
to  one  of  three  different  categories.  These  are:  (a)  routines 
for  proposed  forecasting  or  safety  level  calculations,  (b)  routines 
for  simulating  negative  binomial  requisition  sizes,  and  (c)  modified 
HEDGSIM  routines  required  to  implement  the  above  programs. 

Let  us  now  consider  each  of  these  categories. 


Routines  for  Proposed  Forecasting  and  Safety  Level  Calculations 


This  category  includes  subroutines  FOR576,  LEVELN,  LPGFOR, 
and  SNBROP.  Subroutine  FOR576  is  the  same  as  the  original  HEDGSIM 
routine  with  the  addition  of  an  outliers  test  when  the  forecasting 


code  ICDFOR  *  3.  When  ICDFOR  =  3,  the  subroutine  checks  if  the 
largest  observed  demand  exceeds  the  average  demand  plus  four 
times  the  mean  absolute  deviation  (computed  using  the  remaining 
seven  observations).  If  so,  FOR576  assumes  the  large  demand  is  an 
•'outlier"  that  is  not  expected  to  occur  again  in  future  quarters. 
In  this  case,  the  forecasting  calculations  exclude  the  outlier 
data.  Subroutines  LPGFOR  and  SNBROP  compute  optimum  reorder 
levels  using  the  Logarithmic-Poisson-Gamma  (LPG)  and  scaled  nega¬ 
tive  binomial  models,  respectively.  Subroutine  LEVELN  calls 
these  subroutines  and  also  provides  other  safety  level  computation 
options  within  the  LEVELN  logic.  As  noted  above,  listings  of 
these  programs  may  be  found  in  Appendix  B. 


Routines  for  Negative  Binomial  Requisition  Size  Simulation 

Subroutines  in  this  category  include  DEMPAR,  NEGBIN,  NEGBN1, 
and  RANDEM.  These  routines  implement  the  generation  of  negative 
binomial  requisition  sizes  using  techniques  described  in  Reference  2. 
Subroutine  NEGBIN  utilizes  the  twelve  most  recent  observations 
of  actual  D062  demand  histories  to  estimate  the  parameters  of 
a  negative  binomal  distribution  of  requisition  size.  It  also 
constructs  the  cumulative  distribution  function  (CDF)  of  this 
distribution.  In  turn,  subroutine  NEGBN1  utilizes  this  CDF  to 
determine  the  specific  requisition  size  associated  with  a  randomly 
chosen  percentile  of  the  requisition  size  distribution.  Subroutine 
RANDEM  provides  a  uniform  (0,1)  ramdom  number  stream  which  is 


used  exclusively  in  the  demand  generation  process.  This  guarantees 
that  exactly  the  same  sequence  of  requisitions  is  generated  for 
a  given  random  number  seed  regardless  of  the  management  policy 
being  evaluated.  Finally,  subroutine  DEMPAR  calls  subroutines 
NEGBM1  and  RANDEM  as  needed  to  generate  the  specific  requisition 
streams  associated  with  a  given  demand  history.  Subroutine  NEGBIN 
is  called  to  initialize  the  negative  binomal  requisition  size 
distribution  by  subroutine  INITEM,  while  the  ramdom  number  generator 
RANDEM  is  initialized  by  the  MAIN  program  at  the  beginning  of 
a  simulation  run.  Subroutine  DEMPAR  also  includes  new  logic 
to  insure  that  very  high  activity  items  do  not  cause  the  future 
events  list  to  overflow.  With  the  new  logic,  after  450  requisi¬ 
tions  have  been  generated,  all  remaining  units  of  demand  associated 
with  a  given  quarter  are  placed  on  one  Clast)  requisition.  Only 
extremely  high  activity  items  are  impacted  by  this  change. 


Modified  HEDGSIM  Routines 


Subroutines  included  in  this  category  are  INITEM,  INORD, 
INVRSM,  and  ENTERB.  As  noted  above,  these  routines  contain  very 
slight  modifications  to  the  original  HEDGSIM  programs.  Subroutine 
INITEM  now  contains  a  call  to  subroutine  NEGBIN  to  initialize 
the  negative  binomal  requisition  size  simulation  procedure. 
Subroutine  INORD  has  been  modified  to  generate  GAMMA  distributed 
lead  times  with  a  coefficient  variation  equal  to  .353.  Data 


collected  by  Hayya  (1980)  indicates  that  such  a  distribution 
describes  a  number  of  D062  items  (See  Tables  1-1  and  1-2).  The 
program  INVRSM  is  essentially  identical  to  the  HEDGSIM  MAIN  program 
excepted  it  now  contains  a  call  to  initialize  subroutine  RANDEM. 

When  simulating  high  activity  items  with  average  requisition 
sizes  that  are  close  to  one,  very  large  numbers  of  requisitions 
may  be  generated  within  the  simulation  model.  If  policies  are 
being  simulated  which  offer  very  low  levels  of  support,  this 
can  result  in  very  large  numbers  of  back  orders.  In  turn  this 
may  cause  the  HEDGSIM  backorder  file  to  overflow.  Subroutine 
ENTERB  prevents  this  from  happening  by  cancelling  all  new  back¬ 
orders  once  the  backorder  file  is  filled.  With  the  present  coding, 
once  five  hundred  requisitions  have  been  backordered,  all  subse¬ 
quent  requisitions  are  cancelled. 


Table  1-1 


Coefficients  of  Variation  of  Lead  Times 
Observed  by  Hayya  (L980) 
for  62  High  Activity  Items 


Coefficient  of 
Variation 


Number  of 
Items 


.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 


1.0 


8 

13 

15 

11 

6 

2 

2 

1 

3 

_1 

N  *  62 


Median  *  .36 


Gamma  Probabilities  for 

Mean  =  1  and  Coefficient  of  Variation 

=  .353 

X 

PC  x) 

P(X^x) 

.100 

.000 

.000 

.200 

.001 

.001 

.300 

.007 

.007 

TOT 

.022 

77T30 

.500 

.048 

.077 

.  600 

77577 

715? 

.700 

.101 

.255 

7S7TO 

TTH) 

7T7T 

.900 

.119 

.490 

lTMT) 

Tin 

.602 

1.100 

.098 

.700 

rrm 

7MT 

7781 

1.300 

.054 

.844 

rrm 

75T8 

.892 

1.500 

.035 

.927 

rrm 

77515 

.952 

1.700 

.017 

.969 

T7SM 

775TT 

7987) 

1.900 

.007 

.987 

2.000 

.005 

.992 

2.100 

.003 

.995 

2.200 

.002 

.997 

2.300 

.001 

.998 

2.400 

.001 

.999 

2.500 

.000 

.999 

2.600 

.000 

1.000 

2.700 

.000 

1.000 

Note:  Underlined  values  have  been  incorporated  into  Subrout 
INORD. 
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UP  1 640 , XRS-DEMHY  1HV.0C.H 
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UNITS 

?0,,,10K 
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0 

NOTE 

*.m*«<m>m*«*««IMSM.'ICL  BATA  FILE  FOLLOWS  THIS 
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1 

OPTION 

FQRTRAM,M0«AP 

60 
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SELECT 

INVR/QIJ/1NVRSH.0 

65 

0 

SELECT 

1NVR/0B J/SNBROP .0 

70 

0 

SELECT 

HEBG/OIJ/S INULA. 0 

ao 

9 

SELECT 

HED8/OBJ7REO.O 

90 

t 

SELECT 

HED6/0B J/CUH.O 

100 

i 

SELECT 

MEOG/OIJ/CUMB.O 

no 

* 

SELECT 

MEBG/QBJ/DEHPR3.0 
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SELECT 

HEDG/QBJ/CNTRB3.0 

130 

0 

SELECT 

HED6/0BJ/FILLB0.0 

140 

0 

SELECT 

HEIG/0B J/F  0R573 . 0 

150 

* 

SELECT 

HED6/0B J/IN1 TAL . 0 

160 

t 

SELECT 

HEDG/0BJ/IN1TH3.0 

170 

0 

SELECT 

ME  BG/OB J/ IN0RD3 . 0 

180 

« 

SELECT 

HE06/0I J/LEVLN3 . 0 

100 

0 

SELECT 

HE  BG/0B J/l  ONGSP  •  0 

200 

o 

SELECT 

HEOO/OBJ/OROER.O 

210 

% 

SELECT 

ME8G/0I J/OUTREP  >0 

220 

4 

SELECT 

HED6/0BJ/PFAC.O 

230 

* 

SELECT 

HEDO/OBJ/PLOTR.O 

240 

0 

SELECT 

HEBG/OBJ/RECEIU.0 

250 

0 

SELECT 

HED6/0BJ7RET .0 

260 

% 

SELECT 

HEDG/OBJ/8TATN2.0 

270 

* 

SELECT 

HED670BJ/SSTAT2 • 0 

280 

i 

SELECT 

HEIG/OBJ/ZER02.0 

290 

* 

MOTE 

*«****4****<M*»***t*PREB  routihes  follow - 
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t 

NOTE 
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0 

SELECT 

RE9S/STATUS.0 
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0 

SELECT 

RCQS/FORUPB.O 
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i 

SELECT 

RE8S/LEVEL.0 

340 

i 

SELECT 

REQS/EHTER.O 

350 

0 

SELECT 

REQ8/REH0UE.0 

360 

0 

SELECT 

REQS/URIFEL.O 

370 

0 

SELECT 

REQS/INFEL.O 

380 

0 

SELECT 

REQS/RAN8U.0 

390 

f 

SELECT 

RCQS/6P.0 

400 

0 

EXECUTE 

410 

4 

UNITS 

70,40K, ,10K 

420 

• 

PRHFL 

07/X1SfRf6f 1NUR/I0620C.H 

430 

0 

FILE 

08fX2S 

440 

0 

FILE 

09, ASS 
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* 

NOTE 

**«*»S£L£CTA  NEOG/PFAC.f  FOR  PR06RAN  FACTORS 

460 

0 

BATA 

04 

470 

0 

SELECT A  MEBG/PFAC.I 

INV.OC.H 


480 

• 

DATA 

05 

440 

500 

< 

1  0  t 

510 

0  0  0  0  0 

0  0 

520 

0  0  0  0  0 

OOOOO  0  0  0  0  0 

0  0  0  0 

0 

530 

7  2  38 

540 

2  12  8  1 

2  1 

550 

6  34  0  99 

540 

.20  300. 

500  19500  .5  1 

C7 

570 

1  30  1  500 

8  0  SIMULATE  500 

ITEMS  FOR 

30  UTRS 

580 

2000 

590 

« 

NOTE 

*44*4*4«4»4*4*4*4**4»*«4»4tRINE/G0/PUNCH 

400 

1 

CONVER 

NSPIN 

410 

4 

INPUT 

NNED1A 

420 

1 

OUTPUT 

HODICI 

43  0 

* 

FILE 

IN,A3S 

440 

i 

PUNCH 

OT 

450 

• 

ENIJOB 

INV.OC.H 


20 

« 

IDENT 

UP1 596 , XRS-DE  HMT  INVR.A 

JO 

« 

UNITS 

15,,,10K 

40 

1 

NOTE 

»4«**’MiM*«iM*JNSll.JCL  SATA  FILE  FOLLOUS  THIS 

SO 

1 

OPTION 

FORTRAN, NOMAP 

60 

t 

SELECT 

INVR/OB J/  1MVRSH.0 

65 

0 

SELECT 

1 NOR/OB J/SNBROP . 0 

70 

0 

SELECT 

HEDG/DBJ/S1MULA.0 

80 

1 

SELECT 

MED6/0BJ/REB.0 

90 

1 

SELECT 

HED6/0IJ/CUH.Q 

100 

0 

SELECT 

HE 16/01 J/CUHB.O 

110 

f 

SELECT 

HESG/0BJ/DEMPR3.0 
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t 

SELECT 

HEBG/0SJ/ENTRB3.0 

130 

* 

SELECT 

HED6/0B J/F ILLBO . 0 

140 

t 

SELECT 

HESG/OSJ/FOR573.0 

150 

0 

SELECT 

HED8/OBJ/INITAL.O 

160 

t 

SELECT 

HEDG/0BJ/2N1TN3.0 
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t 

SELECT 

HED6/0BJ/IN0RD3.0 
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0 

SELECT 

HED6/0B J/LEVLH3 . 0 
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« 

SELECT 

HEDG/OIJ/LOH8SP.O 
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t 

SELECT 

HE06 / OB J/ORBER . 0 
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f 

SELECT 

HEOG/OI J/OUTREP . 0 

220 

♦ 

SELECT 

MESG/OBJ/PFAC.O 

230 

0 

SELECT 

HESG/OB J/PLOTR .8 

240 

♦ 

SELECT 

HE06/08 J/RECE 1 V .  0 

250 

1 

SELECT 

HEIG/O8J/RET.0 

260 

1 

SELECT 

HED6/0BJ/STATN2.0 

270 

I 

SELECT 

HEDG/OB J/SST AT 2 . 0 

280 

OfA 

0 

4 

SELECT 

nott 

HE IG/OB J/ZER02 . 0 

Z7V 

300 

f 

i 

nu  ib 
NOTE 

310 

t 

SELECT 

REOS/STATUS.O 

320 

« 

SELECT 

REOS/FORUPD.O 

330 

« 

SELECT 

REOS/LEVEL.Q 

340 

t 

SELECT 

RE8S/ENTER.0 

350 

• 

SELECT 

REOS/RENOVE.O 

360 

t 

SELECT 

REQS/URIFEL.O 

370 

0 

SELECT 

REGS/INFEL.O 

380 

% 

SELECT 

REQS/RANDO'.O 

390 

f 

SELECT 

RE8S/GP.0 

400 

• 

EXECUTE 

410 

« 

LIMITS 

I3,40K,,10K 

420 

t 

PRNFL 

07/X1S,R,S, INVR7I0628N.H 

430 

FILE 

08, X2S 

440 

0 

FILE 

09,635 

450 

0 

NOTE 

♦■♦♦♦♦SELECTA  NEIG/PFAC.D  FOR  PROBRAN  FACTORS 

460 

f 

IOTA 

04 

470 

* 

SELECT A  HEBG/PFAC.I1 

INVR.A 


0  10  0 

0  0  10  1  0  0 

ooooo  00001  ooooo  00000 

7  2  38 

2  12  9  1  2  1 

6  3<S  -99  99 

.20  300.  300  19300  .3  1  C 7 

1  2  1  30  8  0  RUN  2  QTRS,  30  ITEMS 
23 
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SNBROP 
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2  SO  C( )  MMO.N  / 1  O’  )f-()  R/ 1 C  OF  OP 

200  COMMON/  I  CDS  1 0/ 1  CDS  1 0 

2  70  l.W=o 

200  I1'UM=0 

290  I FUM =0 

300  I  RUM=0 

310  KK=NDENf(N) 

320  IF(  KK.GT.  0  )  KK=8 

330  KL*KK- 1 

340  DO  1  1=1, KL 

350  ini)M=I0UM  +  NI)EMNf)(  N  ,  I  ) 

360  IFUM=IFUM  *  NRF0(N,  I) 

3  /0  1  I PUM=I  RUM  +  NRFTUR(N  ,  I  ) 

3  WC 

390  IT=M()!)(  IT  IMF,  I  TOT1?) 

400  R=FLOAT(IT)/FLOAT(ITOTin 

4  10  PM=I  . -R 

4  20  DUM=R*FLOAT<NDHMAC<N  >  )+FLOAT(II)UM  >*RM*FLOAT(NDEMND(  N,  KK  ) ) 

4  30  FlJM*R*FL0  AT  <  N  RF  QAC  <  H ) ) +FL0A  TCI  F  UM  >  ♦  RM*F  LOAT  ( N  RED  ( N ,  KK  > ) 

44o  RUM=R*FLOAT<NRETAC<N)>+FLOAT(lRUM)*RM*FLOAT<NRFTUR(N.  KK  )) 

450C 
4  <OC 
4  70 * 

4HO*i  iCOMPUTF  LEVELS  BASED  ON  NET  DEMANDS,  WHFRF  NET  IS  GROSS 
490**  *  DEMAND  MINUS  SERVICEABLE  RFTURNS  FOR  THE  FAST  24 
500***  MONTHS, 

^  5I0***RETURNS  ARE  ADDITIONS  TO  ON-HAND  ASSFTS. 

FOR573.S 

FOR576 


>4M  IH  I  I  »  A  K..  .  5  »•  ;>  t=  4-i.V  ‘ 

)•>■)  ifo-o  j.i/.o. )  -•)  :  i)  J'hm 

*<*'  f  .U.’Yr,  i  =(  i)  »*/H  . i A  I  (  K<  ) 

.  t )  /i  o  a*j.m  »>  =i  "  «.:•»;  *4 . 

>  v>;; 

>'/O.  J  sjOMFUi1  At  A 1 1  *  >  f-  U>F\  *’!  i  <LY  'i»MAIll> 

;/)>.; 

Jl’l  c.“a>=j. 


*M’5C  ( h  ii  iMvh' 

y-1  1 "  HMC'jI  A'i!)  C\i\\.  Attf-  JA  .*•■)  ON  *1 I  I  OTUl.Y  DIMANU  i<ATf* 

j‘v )  j  : » r* ! •  -,•{ I  iioi'i  i>  A:>HD  on  oi-oss  I ■►-man  I >S 

ofi' j 

‘>/<>  KZI  :!ft)an 

i>JO  DMA  /<=-*>  W)  4  W  . 

•->'/  >  |H)  ?n  1=1  .  KK 

/c> )  DIiM=N  >(  N,  I  ) 

mo  if;(  ii;>m»)c.of'.  j)i>r •  Aar)PM-*Hf»-:ri)N(M. i ) 

/<’ o  IF<  DF  4. OT.  I )  KZM'0=K/H  W)*l 

/*)  I  F(I)F4.I  F.U1AX)  iO  To  lo 
Mi  omaX  =  I)F4 

I  bi )  I  DMA  X=  I 

'OO  lo  CONTINUF 

l  H  •  AUSDF  '/=An.  •  (  I 4-FO  K.\S1‘ ) 

/  hO  CMADeCMA’)  ►AUSDFV 

/  y<r  )  COM  I'll]  f- 

.(■0.; 

no  OF)  -NDFVAOCI  ) 

1^0  I F  (  IIJMOC.';  •.  J)l  F^DFJl-UPPI  JP(N,  I  ) 

1  JO  CMA:)=Ai»  ;(l)*:'4-FO',(;ST*N  )*N 

34  >  CM! )=  ( OMA.'i  +CMA  n )/( 1*1.0  A I  (  K<  )  +R) 

>bO  IFMJMAD.U*.  0.01)  OMAO=O.OI 


no 
‘■<1* ) 

*.  JO 
34  J 
>bO 
i  oOC 
(  /OC 
dHOC 
J^OC 
oyi 
/  IOC 
✓  ZDC 
<  JOC 
MU 
MO 
f  6)C 

^  /oc: 

'  dOC 
WOC 
I  X)OC 
I0IO 

I  o?oc 


IF  i  'OF'CAO'I  COI)F«J,  HACK  OUT  THF  LA WGFST 
UFMANi)  THAT  HXCFFDS  4  MAOS. 

IF<  rCOf'nk.'F-.  j);;o  TO  40 

CIRX  IF  FII.TI-R  IS  TO  IJF  USFI) 

I  F(  KZ':?’O.Li  .  3)00  To  40 
IF  (OMAN.  LI'.  I  .  ).;o  TO  4  0 

IF  DMA  X>YFOl.’  ♦  t*MAO, 

TliFJ  i  SACK  OUT  HIG(?  il  DIMA  NO. 


IFON=<  0.*F  )l?0> T-DMAX  )//. 


FOR573.S 

FOR576 


I 


>’;V 


I030C  CUMPMT'-'  /-IT  MAO 

I  040C 

I  On)  p.  AD=0. 

I  Do  )  1)0  30  1  =  1,  <K 

I  0/0  IK  I  .1-0.  1  -)MAX)  CO  TO  JO 

I  OHO  l'FM-*NI)F.MNO(  N%  I  ) 

I  090  IF  (  II  .  ( IF .  J )  OF  M=|#-  M-  NRf-TtJ R  (N  ,  I  > 

1100  IMAO-TMAD  ♦  AOS(  |>M  -T.  OP ) 

m io  jo  comtpj.jf 
I  120  I'M Ai  l  =  TMA  '/ 

I  I  300 

I  I40C 

I  Iboc  IS  0  UX  AN  OIJI  LI  F«> 

I  I OOC 

1170  I F  ( DAAX.  L::.  <  TFOR+4 .*TMA0  ) ) (30  T()  40 

1  1B0C 

I  19 OC  UbF  OACKFli-O  IT  VALUI-'S 

I  2  OOC 

I P I  0  FOPCST-TFOO 

1220  Al>R(N>  =4.*FoR0ST 

1230  0WAI)=TMAI> 

1240  IMIFHUO.O'-:.  I  )tf  If  II  '*(  6,  33  )I)MAX 

I2SQ  33  FORMAT  (  '  ****F0!?C>  /6--UACK  OUT  OUTLIFR  OF  '.F8.I,'  UNITS') 
I  2oOC 

127.)  40  CONTI  NU,; 

I  2 HOC 

1  29 OC  *  * COMFUTR  STANOAP!)  PFVIATION  OF  LFADTIMF  DFMANI) 

I  300C 

1310  RLT=FI.t)A  r(  LTAI»M<  N  )  +LTR  MOIX  N )  > 

1320  RSIOLTCJ  >=0. S94S*0MAl)*<0. 823/8+0. 4262S  *RLT  ) 

1330  lF(RSh)l.r(N).I.T.  0.01)  RSIGLT(N)=  0.01 

I  34  OC 

I  350C  *  *COMPUTF  AVRi/AGF  PI-OUISITION  SIZE 
I  36  OC 

1  370  IFIFUM.IT.  I  )FUM*I 

I  380  R=*OUM/FiN 

1390  IFIR.LT.  I  .  )R=I  . 

1400  RFi/JIZ(N)  =  P 

1410  CO  TO  2000 

I420C 


/ V.tV.V.VW'm'iVyV.'.'',- 


FOR573.S 

FOR576 


•-UOC  I  ?/»  A  (  .ir-fJI  FOR  Zh IH) -I  AND  H  I  MS 

44  HJ 

4l)  1  |  ■  >00  f.-f; £  ;LT(N» 

•loo  ADMC4  )='>.  ) 

•WO  W  (bi/.(.'()  =  l  . 

4«*)  oo  ro  000  ) 

400  ,O00  IFUFiiUO.IT.  I)  RFITJMO 
boo  WRI TK(Lw.  IOO) 

bl)  WRI  IF  (  LW,  J.  00  )  N,  jMIX:NT( N  ) ,  ICDFO r,  ICDSI 0,  1.  J'A!  M!  N ) ,  I. TP  M0!)(  N  ) 

l>2  t  WRITF!  La,  IOO> 

">  4  »  WRI  TFd.W,  RQo  JNDFMAC!  N)  ,  ( NDFMNIHN.J  )  ,  J=  1  f  KK  ) 

M  >  WRIT1? (  Li.,  2  KO  rUl-l'.u:(  0 >  ,  (  4RH  TO H(N  .  J  > ,  )-  }  t  KK  ) 

bV)  AMI  n:<l  a  ,  ?>  )(>  )NR*-  )A0(  N  > , !  NRFQ!  N  .  J  ) ,  J=  I,  KK) 

Sv>  AMI  TH<  1.4,  4 Of)) 

b  I  )  WRI  TK  (  Lw.ooo)  FORCS  T.  ADR!  M) .  RFOSIZ!  N)  ,  JCOSI  ( N  ) ,  DM  AO,  MSI  GI.T!  N ) 

OHO  100  FORMAT!  I  h  ,  4X  ,  '  **«F()P5 / 3 — N,  NDFNT,  ICIH-OR,  ICDSI  G'  , 

MOA,  '  L  TAOM  ,  LIVM0D7  ) 

oOO  200  FORMAT!  IH  ,  4X,  I  ?  I  //  12  I  /  ) 

olO  300  FORMAT!  Ilf  ,'•><,  'DFM  AND,  RFTURNS,  AND  FRF  O', 

<>?DK  '  It!  OUFMRir  li  IS  TORY  FILF') 

OiO  400  FORMAT!  IH  FOMCST  ADR  AVF-RFO  UCOST', 

04  0K  '  OMAD  M  > I !3.T/  ) 

6b  0  SOO  FORMAT!  I  Ii  ,6FI0.4) 

OOO  PFT1JMN 

o/O  R4D 


FOR573.S 

FOR576 


C.Al'Al.oG/f- It  K  DESCRIPTION^  ilHIO/i 


LNJ.  > 


I  0*«:niN»tHHXj/OilJ/I.cV|.N  1.o(  >CI-,NO';<:) 

2<>*I.FVi-'LN  3.S 

JO  SUHRCkiriNH  LFVm.N(H) 

40C 

50C  TiJIS  ROUTINE  COMPUTES  REORDER,  fiTOCX  OIUI-CTI  VF ,  RET^NTI  ON  , 

o(X:  TFRM  IN  AT  ION  ,  AND  MJ.VoRT  LFVFLS. 

/(x: 

SO  PA  HAM  RTF  i?  NO  JO  =  3J 

ooo 
i  ooo 

I  10  COMMON  / 1  HSUG/IHliUG 

I  20  COMMON/  ICDFOR/ICDFOR 

I  Jo  COMMON/CO  > .IRr/COSHHT 

I  40  COMMON /COS HLU/C03HIJ) 

I  SO  COMMON/CPSOHD/COSOHD  (  J  ) 

I oO  COMMON /Cs  r OPK/CSTBHK 

I  70  COMMON/GoUL  F/G3ULF 

ISO  COMMON  /GM.F/GRLF 

1  90  COMMON /GSLF/GSLF 

200  COMMON/GTL.F/GTLF 

210  COMMON/I riNV/ITINV 

220  COMMON /ITLRVL/n  LEVL  1 

2  30  COMMON/ IOLi-VL/IDLEVL 

240  COMMON /N  ITEM  /N  ITEM 

2S0  COMMON /NDI 1 1 S/NDH  f  S 

2  60  COMMON/POL  IC  Y/I CDEOQ,  I CDSL,  ROOM  AX,  ROOM  IN,  SLMAX.SLM  IN.  RLF.TLF.SUl) 

2/0  COMMON /I  COSl.L/ 1  CD3LI. 

280  COMMON/ ADR/ Al>  R(  I  ) 

2 90  COMMON /LTADM/LTAPM  <  I  ) 

300  COMMON/LTPROtJ/LTPHODd  ) 

310  COMMON /I  ROT  Y/IRQTYOI  ) 

32 O  COMMON/ I RL/l RL  (  U 

330  COMMON/ 1  ROL/ 1 ROL  ( 1  ) 

340  COMMON/ISOI./ISU1.  (  I  ) 

3 SO  COMMON/ISLRYL/ISLFVL( I ) 

360  COMMON/ HI. /I TL(  I  > 

3/0  COMMON/RSI  0LT/R3  IGLF  ( I  ) 

3  SO  COMMON/ RE  OS  I /{/HR  OS  IZ  (  I  ) 

390  COMMON/IJCOS T/IJCOST (  I  ) 

400  COMMON/PF/PF(NOOO), PFA< NQQ0*30> ,ACODF( 30) ♦  NCODR, A I RCRF 

410  COMMON /IISMGC/II SMGC 

420  DATA  Z/I.O/ 

430  DATA  ITLMI N , I ~LM I N /999V9999, 99999999/ 

440C 

4500  COMPUTE  PLANNING  FACTORS 

4C0C 

4/0  AMDA*COSHRf 

4  80  RLTM*LTADM(N  )*LTPROD(N  ) 

400  IF (RLTM.LT. I . )RLTM*O.S 

LEVLN3.S 

LEVELN 


V.'X  V 

^  I ')  Ai  =A  (  H  ) 

*>«:o  6IC  =  KN> 

9  PC 

T‘»fV  l>H  l‘i  ,’M  I  H»-  IM'OGIAM  r -Ai.-J'O'/  i'htHM 

5  hoc 

bf<)  I  OTP-  I J  1  “17 

b  /  ■» 
b»  >C 

V'OO  IF  !<:!H-‘OP  .!■().  V,  A ! )  I !  1 1 ,  |  Pop  POUR' AST  PROGRAM  CHAN? 

6000 

010  If-(  IGDPO'-.  '!!•.  >  >  CO  f  n  >  ! 

G20C 

6J)  PHI  jf);V  =  '  PH  (  I  O  I'  0*  +  fMiUii  ) 

o40G 

OHO  AH  =  PFNOA  *  Ai) 

060  SIG«$IG  *  (  OP NO.V  **  0.85  > 

O  /()C 
60)0 

6^0  20  CONTI NUF 

/OOC 

/ 10  i?LT*AI>Ul.r M/|;\ 

12  o  r/.’p=Ai>/i2. 

/IT)  UOIICDSTOI) 

/40  Af])R=OC*AI> 

/soc 

/ 600  SFT  CO:?  p))IAL  Hi  6MALL  PUHCHASK  OPTM-H  C  )bT 

7/0  COP=COSO’?I><  n 

7  HOC 

/90C  U  RANCH  87  F00  FORMULA  COOP 

bOOC 

HIO  GO  fO  <  I  10,  120,  |  30,  |40,  I  HO,  I  60,  1  /O,  I  80,  I  90)  ,  ICOFOQ 

820C 

8300  2  FIB  /()  C.UIDFLINF  LTD  POLICY 

840C 

850  NO  IF< ADOp.GP.  1000.  )  GOTO  6 2 

860  0=I2.*PMR 

8/0  GO  TO  200 

880  62  IKADDP.GF.5000.  )  50  TO  64 

890  0=  6 . *  PM  P 

900  GO  TO  200 

910  64  0=3.*RMR 

920  GO  TO  200 

9  30C 

940C  AFl.C  S7-o  pul  ICY - DFC  l9/o 

950 C 

96)  120  0*S0'?r(  2.*COP*AI)/(COSHLD*UC  )) 

9  10  AOOST*0*m: 

980  IF < ACOST.L r.CGTO PK )  GO  TO  200 

990  C()R*C()S0WI l)(?) 

I  OO  ■)  Q=SOPr(2.*CORMD/(COSHLD*UC  )) 

ion  go  to  200 

LEVLN3.S 

LEVELN 


020  I  jo  CON  I  I  NU1: 

OiD  140  CON  1  I  NUt;  x 

040  1^0  CON  LI  NUL 

040  I  AO  CONI'  I NUF 

000  I  /O  CON LINUL 

0/0  I  NO  CONTI NUF 

OdO  190  CONTI  NUF: 

090  WRI  TF(  6,8  1 94  )  I CDFOO,  I  CD  31. 

100  8195  FORMAT! '  ****!_  PVR.N —  UN  OFT  I  NI-; )  FORMULA  COOF....', 
II  OK  '  ICI)F.O<}='  ,  13,  '  ICI)5L= I  U 

120  STOP  • 

1  30C 

1 40C  CHFCK  ORDFR  SI/F  LIMITS 

150  200  CONTI NUF 

160  FMX=FO0MAX*RMR 

170  IF(O.GT.HMX)  0=KMX 

180  LMX=!-:OQMIN*HMN 

190  IF(Q.  LT.P4X)  0=FMX 

200  IFIQ.LT.I . )0=  I. 

p  J 

22  OC  B  NAN  CM  BY  SAFFTY  UtVi-L  FORM)  LA  CODF  I  CDSL 

230C 

240  00  TO  <310,320,330,340,350,  360  .  370,380.  390),  ICOSL 

25  OC 

260C  8 Si-'  12  MONTH  SUPPLY 

2  /OC 

280  310  CONTINUE 

290  0=  1 2 . *  RM  P 

300 C 

310  GO  TO  5 00 

320C 

330C  23  AUG  08  OUIDFLINF  LTR  POLICY 

34  OC 

350  320  SL»0. 25*RLT 

360  GO  TO  500 

J7QC 

380C  AFLCM  5/- 6  POLICY 

39  OC 

400  330  Z=SQ  PT (  RHQS IZ(  N  ) ) 

410  GO  TO  358 

420C 

430C  PT-FORMULA  TO  MINI  MI  ZF.  UNITS  BACK')  RDFRFD 

44  OC 

450  340  /.=  ) . 

460  GO  TO  358 

4  700 


LEVLN3 . S 
LEVELN 


I»T— F0UMII1.  A  TO  MINIMI /H  l/»-UU I  >  III u  JS  BACKoR! ) F FFD 
/=RI-.MZ(  N ) 


I  4d  DO 
I  49' k; 

ISO.)  iso 
l  bl  O': 

1  021  )C 

I  o.3()  CONTINUE 

1040  II  (SIO.Lf:.  ).  )0K  =.  ()OOI 

I  000  0010=-!  .411*0/310 

iooo  rnum=ampa*sic. 

I  0/0  IF  ( (0010.  LI.  n.  )  .  AN!).  (0.0 1  ;  '.Ol'.-OO.  )) 

1  080K  Rh  UM  =  A  Mi  >  A  *0  1 0*  (  I  .-  !:x  0  (  OS I  < )  )  ) 

I  09  )  L>N0M=2  .*COOHLI >*UC*/* I  .4  i  4*0 

I  COO  RK=0.  /0/*AL()0(  RNUM/DNOM  ) 

101)  S!.=.RK*SIO 

16?)  00  TO  000 

I  o30C 
I  6400 

I  600C  MODI  F  I  FI)  STD.DFV.  FOR  ACTUAL  OAMMA  U:AD  TIMF  VARIABILITY 

I  oo  JC 

I6/0C  AOSlJMF  COFF  OF  VARIATION  = .  30  3 .  AVI)  I  FFRIOD  »FXPFCTK/)  LFA 

I  OH0C 

10)0  3  00  CONTI  MUF 

I  /OO  FMFANL=I. 

1710  SIOL=. 3S3*RMRAHL 

I  12  0  VAR!  =3101.  **2 

1/30  VARD=b  10**2 

1/40  S  IOOLD=OI  0 

I  /SOC 

I  tone  RFVISF  SIO  FOR  LT  VARIABILITY 

1  /  /OC 

I  /dO  S I G=3G RT (  RMFAfJL*VARD+(  RL  1**2  >*VA  RL> 

!  /9()  IF  ( IMBUG.OT.()).IRITF(  6,-  36  3)SIGOLI> .S  10, S 10  /S IGOU) 

I  done 

1810  3  63  FORMAT  (T20, mS IGOL. D=m , F  10.  I  , M  SIO  LTD  =  '\F  10. 1  ,  ••  RATI0='\FI0 

I  b?oc 

I  83 OC  NOW  JSF  D062  SORT!  /.) FORMULA 

I  84  OC 

IHSO  00  TO  330 
I  860C 

I d 70 ******* ****** ****** **********  ********  ******* 

1  8R  )C 

I89DC  SCALFO  NFGATI VI;  BINOMIAL  MODFL 

I  900C 

1910  370  CONTI NUF 

1  92  OC 


LEVLN3.S 

LEVELN 


I-.V)  I  POP  -  ) 

I  I F (  U).i  V.  )  OO  ro  ni 

I  9‘*0  I P  (  A  /  iA  .  I.  ; . '  i.  )  C.o  i  f)  M  / 

Ivy  K' 

1  V  /■  >C  COmpUTF  CWITIC.AI  ‘.TOOK  Si'OOK  OUT  P  NONA NI  LITY 

I9.M0C 

IVPO  POUT=(  )*C.>  ,il|  l>*UC)  /(  AMI  /A  +a:  >  > 

2onoc 

20  IOC  F  I  :j!  >  (.'(•'[  r  I  CAI.  FILL  PRoiJArf  I  LITY 

?OP() C 

2  030  FCRI  1=1  .-i  UJT 

204  0  I  K  (FCPH  .1  :.C.  J  oo  TO  3/  / 

2  05  DC 

20000  a'O'mgl  nfouisition  size 

20 /oc 

2080  AVI-'-;'-:  J  =  RF‘0:>lZ(N  ) 

2  woe 

?  i  o  oc  eus  T()Mi-i:s  p'-u  period 

2  I  I  oc 

2  120  FD  =  IMR 

2  I  30C 

2I40C  E<LT> 

2  IbOC 

2100  ELT-FI.O AT< I.TAOM <N  >  ♦  LTPROD(N)  ) 

2  I  /OC 

2I80C  ASSvMP  CUFF.  OF  VAR.  OF  LT=.353 

2  IvOC 

2200  CVLT=. 35 3 

22  IOC 

222 OC  FIND  S CAI. FI)  NFG.  BIN.  REORDER  POINT 

2230C 

2240  CALL  SNBCOrM  AVHPEO.El),  ELT.CVLT.  FCRIT .  I  POP ) 

225  OC 

2 26 DC  3AC<  INTO  THE  SAFETY  LEVEL 

22  /OC 

2280  37/  CONTINUE 

2290  SL=I  ROP-RL  T 

2  300C 
231  OC 

2  32  OC  SKIP  SL  LIMITS 

2330C 

2  340  GO  TO  oOO 

2 350*F I LF  FXPFIT.S 

2  360C - 

2  3/OC  EXPONENTIAL  APPROXIMATION  MODFL 

2  3H0C 

2390  3  8 0  CONTINUE 

2400  IF  ( Ai>.  LF.O.  )  GO  TO  38/ 

2410  IMAMDA.LE.O.  )  GOTO  387 

2  42  OC 


LEVLN3.S 

LEVELN 


24T  V  t.'DMi'  l|  !;  Ol'l  lMAI.  STOCKOUT  PNOHA.IILl  iy 

2  44 

24SD  POUT  =  (0*Go,>Hi.|>*UC)/(  AV.DA*AI>) 

24  O0C 

24/OC  l.l.'ii  r  OPTIMUM  HU  I 'Hu  UAH  III  TY 

,'4fiOC 

i-i  ill  i . -H( m.i i 

2SO  XJ 

2SIO  *  I:  t  PI-  I  I.L.  I.  )  r.i)  ID  h  I 

2620  l  H  (  HI- ILL.  ■  L  .  D.UUD'4)  *>!•  I  I  I.  =  (j.  vw 

28  30C 

?S4  >C  AicjiJJIF  !'(/  <=  /H  =  H</.)  =  I  HI  *l;Xl'(  «*/.!> 

2  88  DC! 

2860  /"]  =  u. 

2b /  0  II  C iJl  I  LI  .  D.oc/V/  zr  r.  aL D(»(  (  1  .~H II  I  )/.  331  )/(- 

28U0C 

?t)<K)C  COM2  IT-  M AH,.  AN!)  SAI  l:TY  LFVRL  HI? DM  MAD 

2oO'V: 

261  OC  u  I N’  •  ififOLT  =  .8;48*QmAI>*(.H2  3S  +  )  .42625  *61. TM) 

2  02 OC 

2o  V)  OMAI)  =  Slo/(.-,V45*(.8;M/c+0.42o28*RLTM)  ) 

2  640C 

2  6500  Non  COMPU  i'H  SAH'TY  I  HVRL 

2ooOC 

26/0  SL  *  /r  *  JMAI)  *  b-)HT  (Sl.TM) 

2  66  DC 
26V  DC 

2/OOC  vKIP  THF  I  I M  I  l  CALCULATION 

2  7IOC 

2  /20  CO  TO  600 

2  /  30C 

2/40C  SI I  SI.  ID  2.1-00 

2  /50C 

2/60  38/  CONTI  NUT 

2//0  SL  *  0. 

2/80  GO  TO  600 

2  /yoc 

2800C - 

2810  390  CONTI NUR 

2  82 OC 

2B30C  USH  ’ -PI.  IT-FXPONFNT  I  AL-OAMM A  MODH'L 

2H40C - 

2850  IKAIJ.LT.O.  )  00  TO  39/ 

28o0  If  (AMI)A.LK.O.  )  GO  TO  39 1 

2  8  /OC 

2880C  COMPtJTF  OPTIMAL  SHX^KOUT  PPObAlilLITY 

2890C 

2900  POUT  *  <Q*C06ULn*UC>/(  AMi>A*AP> 

29  IOC 


HOP  2T  s 

463  ) 

WP  HAVR 


LEVLN3.S 

LEVELN 


l.I-'ir  OPilMUM  FILL  PWOMAil  II  I TY 


99200 
2  )  j;  >0 

294  (;  HI- II.L  =  l.-POUT 

2  /woe 

hvoo  n  (Phi  ll. i. !-.o. )  oo  Lo  .19/ 

20/  >  1 1-  (  HP  I  i.l.  :  I' .  0.9999)  PRILL  =  0.9999 

2  9jt  >0 
22-)  )C 

30P0C  FLU  9 A  IF:  OllA WTI-NLY  PHMAIJD  NATH ,  MAI),  AND  LOADIMH 

JO  I  DC 

3 09 00  5INCP  WSlOl.T  =  .  5945*OMA!  '* (  . rt9  35+0  .  426?5*in.TM  ) ,  WF  H  AVF 

30  ioc 

3040  OM A 0  =  SIo/(. 5  >45*(.  3,:  3/5+0.  42625*91. TM)  ) 

3050  n-ATF  -  A  >/4. 

3060  OTWl.V  =  91.  I'M /J. 


30  70C 
30H0C 
J090C 
3100 

3 1  I  OC 
3120 


NOW  OOMPM'IF  WHOWOF  W  LRVFL 
CALL  It  XP  0  A  -i  (  WAT  F,  OM  AD , OTRLT  *PF  ILL,  I  WOP) 
IMIHOH.LT.O)  IWOP=0 


3  I  JO  SI.* I  NOH-Nl.r 

3  I  4  DC 

31500 

JI60C  SKIP  IMF  LIMIT  CALCULATION 

31  /  O'  ‘ 

3I30~  00  TO  600 

3  1900 

3200C  JFT  /OP  TO  ZltHO 

32  IOC 

3220  397  CONTI NUF 

3230  1P0P=0 

3240  SL*I  WOP-OLl' 

3250  00  TO  600 

3260C 

32 /Q*** ****************** ********  ********  ********  ************************ 
3230***** 

3290  500  CONTI  NUT: 

3300C  ! 

331  OC  BRANCH  BY  SAFFTY  LFVRL  LIMIT  CODR,  ICPSLL 

3  32  OC 

3330  00  TO  (‘>10,520),  ICDSLL 

3340C 


3350C  D0<>2  LIMITS,  AS  OF  JUNR  // 


3360C 

33/0  510  CONflNUF 

3J30C 


LEVLN3 .  S 
LEVELN 


<33  >r 

3-iooc 

34 1  cn: 

LIMIT  ‘«AFFTY 

i.,;V  l:L 

)C 

3  1 30C 

LIMITS  »!.  TO 

LF1'^  THAN  L»;A01TMH  Dl-MANI) 

3440 

5LM  =  N|.  1 

3450 

3460C 

IF(SL.«M  .  >  1 . M  )'iL  =s 

LM 

34 /OC 

3  A  HOC 

1.13!  1  SI.  TO 

LFSS  1  HAN  3*0 !  •  ’> 

3490 

LLM=3.  *SI  I 

3500 

35 1 0** 
3520C 

IF(SL.<;r.  >1  M  JS1.=S|.M 

35  30C 

3540C 

1.1  MI  !S  UASR) 

ON  MO  NT  it*  OF  SUPPLY 

3550  5  20 

CONI  1  NIJI- 

35oOC 

i:oN'-.<  LIMIT 

35/0 

3LM=8I.M|  J*KMW 

3580 

IFlSL.Li V  LM  ) S 1  - 

SI.M 

3590C 

JPPFR  LIMIT 

3000 

.8  |_M =S  LM  A  X  *  RM  R 

3610 

3  62 OC 

3630C 

IF(SL.OT.SLMKSL= 

81.  M 

3640 

GO  TO  600 

3 6bOC **************  ********  ******  ******************* 

3660C 

30/0  00 0  CONTI NUF 

3O80C 

3690C  COMPUTE:  Lh'VFLS 

3700C 

3/10  ISLf-:VL(0)=SL+0.*> 

3/20  I  RUTY (  N )  =0+0.5 

3/30  IR0L<rn=5L+RLT+0.5 

3/40  IFCIPOL(N).LT.O)  IROL<N)  =  0 

3750  ITL(N)  =S  I.  M  A  X  *R  MR  +1?  LT  +OTLF*  RMR  ♦  0 . 5 

J/oO  1 F  ( ITL  ( N  )  .  l.T.O  )  I  ]"L  (N  >=0 

3  HO  I  RUN  > -FLOAT <  lTL<N>)+GRLF*RMR+0.5 

3/80  IF  <  ITL  <  N  > .  I  F.ITLMIN  )  IT1.C  N  J-ITLMI  N 

3/90  I F  ( I  RL  ( N )  .LF.I  RLM  I N  )  I  RL  ( N )  *  I  RLM  I N 

3800  ISUL(  'I  )=GSULF*RMR+0.5 

3810  3000  CONTI NUF 

3820  ROL -FLOAT ( I ROL ( N  ) ) 

3830  IF( IHBUG.LT. 1 )  RFTURN 

3840  RK  =0. 

3850  IF(SIG.OT.O.  )  RK  *  SL/lilG 

3860  h'RI  TF(  6,8903)N*  IROTY(N>  ,1  ROL(N)  ,  ITL(N),  l RL( N) , I SUL< N>  ,  RK, 

3  870&  SL ,  PFNOW, AD.PF I  LL 

3880  8903  FORM  AT  <  4  X,y  ***★!.  FVFLN--  N*'  f  15,  '  I  ROTY='  ,  15,  *  IR0L=y,I5, 
3890A  '  ITL=',I8,'  IRL= 15, '  ISUL='.I5,'  RK=',F9.4, 

3900A  /,T50, '  liL=',F9.1f'  PFNOW*'  ,F6. 3,  *  AD-y  ,F9.  I 

39 1  OK  PFILL  -  ' ,  F  6 .  3  ) 

3920C 

3930  RFTURN 

3940  FND 

LEVLN3.S 

LEVELN 


SUMtttiUriNI  r;XPL  1 IX  <.  W,QMAI),OTrtl.T ,CUMPX  > 


4  /*-,■ , 

4  /f> .  >9 

•i//  ;  common /|Mr/IftT( 20) 

4  /d  ) C 

4 /V  JC  SET  WRITE  FI  AOS 

4  HO  >C 

4  HI  0  II'FTL  *  K.T(IO) 

4H20  IPNTS/  «  I WT( I  I  ) 

4«d()  CUMPT=0. 

4840  CliMPX»0. 

4  850C 

4 9600  INITIALIZE  PDF  PARAMETERS 

4870C 

4883  A I  =  ').  3  II 

4890  Bl  *  -0.4H.3 

4  9000 

491 !)  A 2=  0.  009 

492  0  82  =  0.  /9/9 

49.300 

4940C  GAMMA  CONSTANT  FOR  MEAN*  I  AND  CV=.  3 53 

4950  Cl=  0.  OOI5H/3 

4960C 

4970C 

4 9800 - 

4990C 

bOOOC 

501 OC  INITIALIZE  FOP  T  INTEGRATION 

50200 

5030  l)T  *  .1 

5040  T  =  DT 

5050  CUMPT=0. 

5060  CUMPX=0. 

50 /OC 

5080C  BEGIN  "TM  INTEGRATION  L(X)P 

509 OC 

5100  DO  100  1  =  1  ,  I  00 

5  I  IOC 

5I20C  COMPOTE  STANDARDIZED  ERROR  ZT 

5I30C 

5140  TCTR=  T*ar?LT 

5150  ZT  »  (X  -  R*TQTR>/  ( QMAI>*SQRT<TGTR )  ) 

5  I  60C 


LEVLN3 . S 
EXPLTD 


b  I  /  ■ )  ‘  I  l;  r*  (  I  ) 

b  i  > )  cr;  =  ci  **/*  •  <p(-m.*t>  *  dj' 

b  2  />  »o 

•jjHc  comp.i  1 1-  p<  /.  <-;  i  \  n 

b  .■> 1  >e 

->23o  IK(  /.I'.l.K  )  P'  -A  >*t  Xt'(ll2*/.T) 

ik  /r..i;r,  j.)  i*/.=  i.  -  ai*kxk  iu*/.r> 

*■)  ">■>  ■  )C 

■>?c  )C  C OMI’I.’TK  h  X  X  11)  K  i) 

b2/oe 

b?n  )  i>xt  =  p/*  ; . 

“.>2  9  )C 

b  JO  )C  tJPOA  n;  CsJMUl.Ar  I VM  PKOUAUl  L  ITIFS 

b  JlnC 

bjp'>  cuMPi"  =  cump  r  +  or 

b  33 )  CUMPX  =  CUWP.X  ♦  HXT 

b  34  y  : 

b  3b  O  IPPNT=0 

b  Jo  )  IK  (IPNri'.CI  .0  )  .AND.  ( MOD  (  I .  I  PNTSZ )  .  f-Q.  0)  ) 

bj/)x  iPw,Nr=i 

b  380  IMMKTL.LK.O)  I  PPNTO 

bJ‘>o  IK(IPPNT.Ol.O)  KvPin-(6.63)  X  ,T,  Gl,  PXT, CUMPI, CUMPX  t  ZT,P£ 

b^fn  oj  R)PMVr(2F  i.?t<')K  J0.4  > 

o4i  k; 
b  4?')C 

b4300  I NOKFMFNT  T 

b  44  OC 

b4bO  T  *  T+DT 

b  4oOC 

b4/0C  STOP  IK  CUMPI  >  .999 

b4HOC 

b 490  IK (CUMPT.OT.  0.999  )  00  TO  1?0 

bbOOC 

bblOC - Rll)  OK  "T"  L(X)P - 

5520  100  CONIINUK 

bbJOC 

bb49  120  CONTI N UK 
b550  Pf-TUWN 

bb60  FNO 


LEVLM3.S 

EXPLTD 


W50  SUiJKOUriM!-  FXP  JAM(  OPAl'F,  WAI).  TOTR.PL  IL.I  .  I  QoP  > 

j  ;ooc 

n/OG  COMPU  l'i;  A  HF-OHDF.I?  HO]  Mi  (If, ’Op)  T:)  0 1 V I- 

3  >8 >0  a  f:  III  PAIR  of  PL  III.. 

OC 

4  000  COMMON/ IWIV  I  WI  ( 20) 

40 10(; 

402X'  SFi  LOW  HO  I N  T  FOt?  P.I).  L 

4  0J0C 

4040  QMAOLT  =  :)MAI)*bOPf  <  TOTR) 

4050  Xl.OW  =  QP4  TF*ToTR  -  5.*QMAI)LT 

4060  XLOW  =  IITXfXLOP  +  0.5) 

40/0  I F  (Xl.OW  .LT.O.  )  XLOW  =  O. 

4 on  CALL  FXPLT  )  (  XLOW  ,  ORATF.OMAI),  TOTP.CIJMPX  ) 

40V0C 

4  100  IFfPFILL.  ii'.ciiMPX)  (JO  To  500 

4  I  I  OC 

4I20C  VFPY  LOW  SUPPORT  HFFOFI). 

4  I  3  )C  S FT  '(FOIL M  R  POINT  TO  XLOW 

4  1 40C 

4  150  I  POP  =  XLOW 

4160  RLTU7N 

4  1  /  OC 

4I80C  SFT  LOW  POINT  FOP  SFARCH 

4  I90C 

4200  500  CONTINUF 

4210  PL()W=  CUMPX 

4  220C 

4 230C  COMP JTF  OFLTA-X 

4240C 

4250  OX  *  GMADLT 

4260  II  (OX.LT.  I  .  )  l)X  =  I  . 

4270  IF(DX.OT.I  .  )  l)X  «  IF IX <  l)X  +  0.5) 

4  280C 

4  290C  INITILTZF  HIGH  VAI.UFS  FOR  SFARCH 

4  300C 

4  310  X*  XLOW  ♦  l)X 

4  32  OC 

4  330C  FINO  PCDI.T  <=  X) 

4  340C 

4  350  510  CONTINIJF 


LEVLN3.S 

EXPGAM 


4  30)  CALL  FXFl  .1'  X  X.  -WAT'-,  -ml’ .TOTI  V'UMPX  ) 

4  3/1  IF  ( IV<T(  IO>.OT.-i  )  W  MT';  ('>. 523  )X,C!)MPX,X  LOW.  PLOW 

43«)  ‘93  FORMATC" - EXPOAW.  X*",Fd.l,"  CUMPX=", Ed.4 , 

4  39 )*.  ••  XL>  >;Y  =‘*  ,F>i.  I  HOW  ="  ,FH.  4 ) 

44(,’)C 
•1  4  I  >C 

44  9)C  HA/1'  WE  MOUNDED  flH-  IX-SI  RED  Fill.  WAIF 

4430C  IF  so,  :;o  To  oQO  AN  I  >  IN  T  FRPOLATF  FOR  R()F>. 

444  0C 

445<)  I E ( PE  I LL  »LE.  CdMi'X  )  <>0  TO  oOO 

4  4O0C 

44/OC  IF  OJVIPX  >  .99,  bToi*  ANYWAY. 

449  )C 

4490  IF  (CJMPX  .0'-'.  0.99)  00  TO  590 

4500C 

4 5 I OC  NO,  <E  HAVF  TO  KEEP  TRYING.  NF‘>F T  FOR  SEARCH 

4520C 

4530  XI.OW  =  X 

4540  FLOW  =  CUMFX 

4550C 

456  )  X=  X  ♦  1)X 

45/0  00  TO  510 

4  56  OC 

459  >C  SFT  IROP  =  THIS  X 

4O00C 

4610  5 HO  CONTINUE 

4  620  IROP  =  X*  0.5 

4  030  RETURN 

4640C 

4  650C  THE  ROP  IS  HOUNDED.  NOW  INTERPOLATE. 

4600C 

46/0  600  CONTINUE 

4  69  OC 

4690  X  =  XI.OW  ♦  CPF* 1 1.1.  -  PLOW)*  (  X  -  XLOW)/  (CUMPX  -  PLOW) 

4  /OOC 

4710  IROP  =  IF  I  X(  X  ♦  0.5) 

4/20C 

4/30  RETURN 

4  /40  END 


LEVLN3 . S 
EXPGAM 


CAl'ALOO/F  ILF  DESCRIPTION 


Invr/lpgfur.., 


s*fim  ipgfor.r 

0  GC 

tc  FORTRAN  V FUSION  of-  LPOFOR.b 

0  8C 

I  oc 

IS  CHARACTER  ADol.R*  I 
I  6C 
POC 
PPC 

25  CHARACTE R  Ai)f)LR*I 

30  COMMON/ LPGPRM/T  1  ,  RL  I  *  A I  .til,  I  DM  1)8 

40C 

oOC  T(  K+l  )  =  T(X,K)  FOR  THE  CURRENT  X 

OOC  TP( K+ I  )  =  T(X-I.K)  F RO M  PREVIOUS  X  CALCULATION 

/OC 
HOC 

/OC  REAIH5,  I  3)  P  ARAMFTFRS 

9  1C 

9 2  90  CONTINJF 

9?  PRINT  ,  "IDY.IDH  •' 

93  Pf-AD(5,I3>  1 08,  ID/ 

V  9  AC 

98  13  FORMAT  (V) 

100  CALL  LPGO 

I  IOC  00  l.PG  RECURSION 

120  CALL  LPGREC 

130C  DO  SCALED  NEG  BINOMIAL  CALCULATIONS 

140  CALL  LPG3 

I  42  C 

I 43C  CALL  NAHMIAS  EXACT  LPG  ROUTINE 

I  44  C 

145  CALL  LPG  I 

I46C 

150  WRITE!  6,  13)  "CONTINUE  (Y  OR  N)M 

155  REA 0(5,  I  3)  ADOLR 

I  00  I  F(  ADOLR.  EO.  "Y"  )G0  TO  90 
I  /O  STOP 

I  71  END 

I  /PC 

I  /3C 

174  SUBROUTINE  LPGO 

I  75  COMMON/LPGPRM/TI ,  RL  1,  Ai  , B  1 ,  1 07,  1 1)8 

177  13  FORMA T(  V) 

I  /8C 

1/9  DIMENSION  T ( 200 ) , T2  ( 200  ) 

I  HOC 

181  WRITE(6,  13) 

182  WRITE  (6,13) 

190  WRI TF( 6,  I  3)  "LPGO.S  — FXACT  LPG  PROB  CALCULATIONS  USING  RECURSION4' 

191  W  RI TE  ( 6 ,  13) 

192  WRITE(6,  13) 


LPGFOR.S 


ai 


N  •-  -- 


200  ' 

210  l-‘,I  i  K(  o,  1  3)  ‘"“ADC),  IJ>  AVI-  RFC,  !-(!>>,  H/U>,  C.OF.V  OH  L  T 

?>0  i,FAP<5,M)  P) . P 1 ,  f  I  ,0 
2  30  WPI  TF  (  3,1  3>  PO,  1)1  , H  I  ,c 
2400 

250C  b(^L'/l:  F04  TiJH'l'A  =  Tl 

2  CO  CALL  JS'fCH  (  RO,  0 ) 

2/0  T  1  -Q 

280C 

2900  ESTIMATE  PARAMETERS  FOP  I. HO 

3000  S  =  STD  PH'/ ,  13  =  VAP  TO  Ml- AN  PATIO 

310  S  =  C*!:l 

320  3  =  S*3/E I 

330  U  1  =  l)l*FI 

3400 

350  131  =  1/3 

3  60  AI  =BI*^I 

3/0  R  I  =-T  I  //  <  I  -  T I  )  *AL<X>Z  I  — T I  i  ) 

31 30  RLI  =  Ul  /<  EI*PI  ) 

390C 

400  W  R I  TH  ( 6 ,  I  3  )  •*  A  V F  RHQ  3 1  ZF  =  ",  PI 

4  10  Wf’I  TH C  6,  I  3) 

420  WR ITF  /  c>,  |  3  )  "THETA  =•* .T 1  ."LAMBDA  =  *\RL1 

430  WRITE  /  6,  I  3)  "  ALPHA®"  ,  A I  ,  MiJHTA  =  ",B1 

440C 
450C 

460  RMI=-TI/<(l-n  >*ALOGI  I -T I ) *RL1 *AI /Bl  ) 

4  /O  WPITE/6,  I  3) 

480  W  PI  TE  (  6,  1  3)  "MEAN  =" ,  RM  I 
490  WRITE/6,  I  3) 


500C 
510 
520 
530C 
540C 
550 
5  60 C 
5/0 
5800 
590 
6000 
610 
620 
630 
6400 

650 

660 

670 

6800 

6900 

too 

no 

720 

730 

740 

7500 

7600 

.7700 

780 

790 

800 

801 


ESTIMATE  THE  FIRST  FOUR  MOMENTS  RM I ,  RM2,  RM3,  RM4  OF  THELPG  DISTRI 
d i  *  ui*/  1-ri ) 

04  *  -PL  1 *TI // ALOGZ I ,-TI ) ) 


A  I  *04/1)1 


RM2  =  A  I  *04*/  B I  +04  )  /D 1  **2 

RM3  =  A 1  *04*  ( / 131  **2 )*/  I  *T  I  )  +2*C4**2 )/D 1 **3 

RM4  *  13 1  **3*  /  1  *4*T  I  *T  1  **2 )  +  BJ**2*C4*(  3.*AI  «■  1.) 
RM4  =  RM4  ♦  6*131  *C4*C4*AI  +C4**3*  (  3*A  I  *6.  ) 

RM4  =  A  1 *C4*RM4/DI **4 


RM4" 


WRITF/6,  13)  -"MOMENTS* 

WRITE/ 6,  l  3)  "  RM1  RM2  RM3 

WRITE/6,13)  RMI , RM2, RM3,RM4 

COMPUTE  STANDARDIZED  MOMENTS 
WRITE/6,  13) 

S  =  SORT/  RM2) 

WRITE/6,13)  "COEF  OF  VAR  S/PM  1 
WRITE/6,13)  »RM3/S**3«  ",  RM3/S**3 
WRITE/ 6, 13)  «RM4/S**4»  ",  RM4/S**4 


OUTPUT  LPG  PARAMETERS  TO  FILE 

WRITE  18,13)  RMI  ,S/RMI  ,RM3/S**3,  RM4/S**4 

WRITE  / 8 , 1 3)  Tl ,  RL1 , A  1 ,B I 

RETURN 


LPGFOR.S 


809 
H  MC 
M  >.)C 
8  3  *C 
d  JIC 
833 
8  El 
6  40  ; 

85')C 
80')  C 
8  O-j 

8  /o 

8fX) 

890C 

yooc 

y  ioc 

920 

9  30C 
9  40 
V50C 
960C 
9  70 C 
980 
990 

1000  H I =S7*H I 

IOIOC  S/  *  SCALE  FACTOR 

1020  H2  =  (C/!  RH+Bl  >> 

I  030C 

1040  W  HI  TE  <  6 , 1  3  >  "C ,  C  1 ,  H 1  * H24\C  f C  1 ,  H 1  , H2 
104  1  13  FORM  AT  C  V  ) 

I  05  OC 

I060C  evaluate  H!X=0) 

10 /OC 

I  080  Til)  =  O 
I  090  T2( I )  =  HI 

MOO  51  =  ,11 

1110  C2  *  C2  ♦  SI/S/ 

I  1 20C 

1130  WRITE (6,  13) 

I  131  WRITE!  8,13)  0,SI  ,C2 

I  I40C 

1141  WRITE! 6, 13) 

1150  WRITE!  6,  1 3)  ••  X  H!X>  F( X )" 

1151  WRITE!  6,  I  3) 

I  160  WRITE! o, I  3 )  O.SI ,  C2 
I  I  /OC 
I  1 80C 

I I90C  EVALUATE  H!X)  FOR  X  .GT.  0 

1200  1200  CONTINUE 

1201  L2»RL2 

1202  00  1050  IX= I , L2 

1203  X*IX 

1204  IXPJ=IX-H 

1220  SI  =  O 

1230  T(J)«0  LPGFOR.S 

-  LPGREC 


s  iti  ’our  i  oi  i.i'f  i  i-v 

LEG  RECURSION  C  U.CULA  1  I  ON 

★  A-*-****  V*  •**  *  *  *•*  ♦****)■  *-*•*-****★★•**★**★*  ★★  kit  kit  it  kifkirirk 

COMMO  1/LP3C  >‘f/[  I  ,1)1.!  ,AI  ,3  I  ,[:>/•  I  DM 
1)1  Mf'NSlON  1  (  MO),  I  ?!  >00) 

’EUEFiNI-  C  FOR  RRoTAOIi  nV  CALCULATIONS 

WRITE!  o,  I  D  M  A*A-l.Rr.  .;i  c,  H  ,  l,L  I  ,AI  ,  81  ".Tl. PLI.AI  .81 
C=  -  )l.l  /Aid'd  I  .  —  II  ) 

C I  =  C/<  ELI  +  81  ) 

SET  Li  MU'  ON  X  =  2  00 

31.2=20  0 
C?=0 

SET  CONSTANTS  FOR  USE  IN  RECURSION 

S  /=  1 

H I = ! B I / ! RLI+81 ) ) **A I 


✓4.x; 

25)  it  1  X H  )  =  ITI*H2*(AI  +X-I  )/X)*T2(IX> 

2000 

2/0  IK  X  .LI.  2  ><50  TO  I  OHO 
2nOC 
2  V  )C 
300C 

310  IX)  I  16')  K  =  !  ,  1 X—  I 

020  T(K+)  )=tTI/X)  *  ( H  2  *  ( A  I  +  K-  I  )  *T?  (  K)  +  (  X-  I  )  *T?(  K-H  ) ) 

OJO  SI  =SI +T(K+ I  ) 

340  I F (  10  i  .LI-..  0  >00  TO  I  3 60 

3b  0  WMl'Eto,  IJj“X,K,T(  X.IO.SI  ,  X ,  < ,  T  (  K  ) ,  S  I 

36  .)  l  J60  CONTI NUF 

3/3C  PICK  UP  KX.X)  TERM  IN  SUM 

3/S  I  3 HO  CONTI NUK 
JrtO  SI=SI  ♦Ttlt+IK 

39  OC 

400C  WRIi'h'tO,  13)  TOTALS  FOR  H (  X ) 

4  IOC 

420  C2  =  C  2  +  SI/S/ 

430  WRITE(6,|J)  I  X.  S.l  ,  C2 

440  WRITE  (8,13)  IX  ,  S.l  ,C2 

4S0C 

4O0C  I  Ft  CUM  PROS  FXCf-FDS  .99,  STOP 

4/OC 

4 HO  IFt  C2  .GT.  .99  )GO  TO  1670 

490C 

SOOC  WRITE (6,13)  TtK)  TERMS  FOR  DEBUGGING 

SIOC 

S20  IFt  1 08  .LF.  0  >00  TO  15  60 

530  LX)  1550  KK=I,IXP| 

535  K=KK- 1 

540  WRITEt  6,13)  "X,K,T(X,K)  »>M ,X,K,T(  KK  > 

55 0  1550  CONTINUE 
560  15  60  CONTINUE 
5/OC 

580C  RECORD  T(K>  VALUES  FOR  USE  IN  NEXT  PASS 

59  OC 

600  DO  1620  KK=1 ,IXP| 

610  T2(KK)=T(KK) 

620  1620  CONTINUE 
63  OC 

64 OC  - END  OF  X  LOOP 

650  1650  CONTINUE 
660C 

665  16/0  CONTINUE 

670  WRITE  (8,13)  -99,-99  ,-99  END  OF  LPG " 

680  RETURN 

683  END 

&?0C  *********************** *★*★***★*★*★ ****FND  OF  LPG  RECURSION 


LPGFOR . S 
LPGREC 


FVT.'T 7 


,  *  »  ■  .  X  7* 


i 


i 


•v,  ' 
%  % 


i: 


p. 


N 


/. ; ) : 

/i  ' 

/:*  k: 

i  .>  i 

/  i  ij 

/  4  •  )C  ■ 

74 -1C 
/4)C 
/4fXJ 
/‘ioc 
VfOC 
//  0  c 
/8  )C 
/V  JC 

dcnc 

hi  k: 

h2  )C 
830C 
«40C 
85  O 
880 
H/0 
HHO 
890 
I  900 
9  IOC 
I  92 OC 
9  iOC 

>35  1940  CONTINUE 

94  0  I  F  (  RO  .0  ; 

95  oc 
9  CO  C 

19/0  R9=H 

980  Q9=Q 

1990  GOTO  2080 

199  1C 

992  2000  CON  r  I  NUH 
2  OOOC 
2010c 

2020  RI=R 

2030  01  ='* 

2040  GOTO  203) 

2  05  OC 

20o()C  HALVE  THE  INTERVAL,  ANO  RE-EVALUATE  THE  FUNCTION 

20 /OC 

20/5  2080  CONTINUE 
2080  K=K+I 

2090  0= (01+09 )/? 

2  100  R=FNR(  0) 

2  I2  0C 

^|30C  IF(  R  Id  WITHIN  ,001  OF  Ro,  )G0  To  RETURN 

2  I40C 

2141  IF  (  1 1)8. GT.  ) )  WRITE ( 0,2  3/  K,RO,R,Q 

2142  23  FORMAT  (  "  ***HSRCH~K,  RO  ,  R,Q‘M5 , 3F8.  3  > 

2150  IF(  AHS(R-RO)  .LT,  .001  )G(>  TO  2180 

2loO  IF(  K  ,CT.  15  )G0  TO  2180 

21/0  GOTO  1940 

2175  2180  CONTINUE 

2180  RETURN  LPGFOR.s 

BSRCH 


)u  T  I fi- -  A:,  \;n(  R  ),  ) ) 

COMMO  l/u>  JIH/M/n.lJ  I.M.h:  .  !'  /,  III.) 

- -  v(.(j  I  /!-;  fouAllOM 

-  KJN.TION  FNR(f»  =  --  /(  (  I  -  0 >  *  Al  O  !  (  |  - O)  ) 

•'FF  >102-2190  F >>R  FN7(f.O  )  I  I  I  5'  IT l ON. 


HI  NARY  M-A'-C  I  ROUTINE 


301. VR  FOi.’  THE  VALU'-  0  WHICH  G|V1 \>  AND 
AVI-'.  RFCIJISI  n  ON  }f/|-  Of-  i’o 

SET  UR  END  POINTS  FOR  :>FARCH 

K  =  O 
09*. 99  > 
i?9«FNR(  09  ) 

01  =.001 
Rl  =FNR(QI  ) 

GOTO  2080 

CHECk  I F(  RO  .LF.  R 


\  R  )GO  TO  2000 
RE-SET  TOP  OF  INTERVAL 


RF-SFT  80TT0M  OF  INTERVAL 


2181 


END 


2182  FUNCTION  FMIHO) 

2  183  FNP=-Q/<<  l-Q>*ALO)(  l-O)  > 

2184  PFTUPU 

2185  FND 


LPGFOR.S 

BSRCH 


2  380  SUBROUTINE  LHtJl 
2400  I  <1  MANSION  Y  (I  00,100) 

2410  COMMON/LPGPPM/l'l  ,  PL  I  ,A  1 ,  ill  ,  in/,  11)8 

2411  WRITF(6'M>  “T  UG  PROGRAM  COMPUTES  EXACT  PIft)M  FOR  IMF  IPG  GIST" 

2412  »3  FORMAT!  V) 

2420  W7IfF(6,l*>  "  TH  FT  A  =  M,T1  ,  "LAMBDA  ="  ,  PI.  I  ,  "ALPHA  =  ",  Al  ,  "BFTA  =  J,,BI 
2430  WRITE! 6,  I  A)  I  ,2, 3,4 

2440  WRITE!  o,  I  3>  I  I ,  IM.t  ,  A  1 , 31 

2450  C  «  -  RU  /Al  00(  I  -Tl  ) 

2  460  PM  I  =  -TI/((I-I'I  >*ALOG(  l-TI  )  )*PLI*M/8I 
24/0  Cl  =C/(  EL I  ♦  i)|  ) 

2430  L  2  *  I  FIX!  I  OO*  EM  I  ) 

2490  C2  =0 

2500  I)!)  2/90  I  X  8  I  3 1  ,*l.  2+ 1 

2504  X=I XP I- I 

250  6  I  X=X 

2510  Y(  I  XP  I  ,  I  )=0 

2520  Fl-I. 

252  1C 

2522C  COMPUTE  !X-I) 

252  3C 

2530  IF!  X.LT.2.JG!)  To  25 /O 

2540  DO  25oO  K=*l  .  IX- 1 

2550  FI  =  F I *K 

2560  2560  CONTINUE 

2565  2570  CONTINUE 

25/0  Y  <  I  XP  I ,  I  XP  I  )«I./FI 

2580  FI  -  I 

2590  IF!  X.GT.  O.  )Oo  T!J  26  30 
2600  SI  =  !BI/!ELI  ♦8I))**AI 
2  610  GO  TO  2740 
2620C 

2625  2630  CONTINUE 
2630  Sl=0. 

2640  DO  2/20  KK=I,IX 

264  1  KP I  =KK+  I 

2650  IF!  IX  . EQ.  I  )G!)  TO  26/0 

2660  Y!  IXPI  ,KPI  )=Y!  I  XP  I  —  I  ,  KP  I  —  I  ■)/!  I  XP  I  —I  )  +  Y  !  I  XP  1  - 1 ,  KP  I  > 

2665  2670  CONTINUE 
26/0  FI  =1  . 

2  680  DO  2700  .JJ  =  ! ,  KfC 

2085  J*J3-I 

2690  F  I  *  F I  ★!  A 1  ♦  J)*CI 

2700  2/00  CONTINUE 

2/10  SI  =  SI  +  YIIXPI  ,KPI  )*F  I 

2/20  2720  CONTINUE 

2 730  SI  *  S  !*( B I  /!  RL  I  ♦  B I  ) ) **A l*TI  **X/X 
2/35  2/40  CONTINUES 


LPGFOR.S 


274  )  (  2  =  C2  +  SI 
2  A'O  IF  (  ID8  .I.E.O  )()()  TO  ?/80 
2/60  W.'iITE(6,I3)  X ,  i>  I  ,C2 
2  HO  IM?  ITF  <8. I J>  X ,  S  I  ,C? 

2 //I  2/80  CONTINUE 

2/80  I F  (  C2  .OT.  .9/  >00  TO  2800 

2  /VO  2790  CONTI NUF 

2798  2800  CONTINUE 

2800  W.'/ITE  (8.13)  -99.-99,-99 

2810  PE  TURN 

2815  END 


LPGFOR.S 

BSRCH 


2h2')C 

2  830C  ★***  ****  A-*********  *★★*  **  *****  *  ********  **  **********  ***** 

2840  SUBROUTINE  LP  )2 

2s4s  COMMON /LPGPRM/T I  ,RI.I  ,  A  1 ,  lit  ,1!)/,  ]  08 

2850  WRITE(6,  1.0  “SCALP)  POISSON" 

285s  I  3  FORMAT! V ) 

2  8oOC 

28/0  WRITF(  6,  1  3)  "1'IFTA  =  '\T1  ,  "LAMBDA  =" ,  RL I  .  "ALPHA  =  ",AI  ,  "BETA  =4',BI 
2880  C  =  -RL1/ALOC.(  I- II) 

2890  C 3-  I./( I • -Tl > 

2900  K I  =  T I *C 

2910  C4  =  <BI/(KI+3I  )  )**AI 

2920  A  =  KI/(KI+BI> 

2930  C2  =  0 

2940  00  3090  MM=I,IOOI 

2945  M=MM- 1 

2950  PI  =  I 

2960  IF  (  M.NF.O  )G0  TO  3000 

29/0  PI  *  C4 

2980  GOTO  3040 

2990C 

2995  3000  CONTINUE 
3000  IX)  3020  11  =  1  ,  M 

3005  I=II-I 

3010  PI  =  P I  * <  <  A  I  ♦M-I-l  )/(M-I  )*!V) 

3020  3020  CONTINUE 

3030  PI  «  PI*C4 

3035  3040  CONTINUE 

3040  C2  =  C2  ♦  PI 

3050  IF <  ID8.LE.0  )G0  TO  3080 

3060  WPITE(6, I  3)  M, C3*M,  P I f C2 

3070  WRITE  (8,13)  C)*M,PI,C2 

30/5  3080  CONTINUE 

3080  I  F(  C2.GT..99  )0()  TO  3100 

3090  3090  CONTINUF 

3095  3KX)  CONTINUE 

3100  RETURN 

3101  END 


LPGPOR.s 


★  ***•  4r***  *•*  ****  Irk  k  k  ***  *  kkhkk  »  *  •*  k  k  kk*k k  k  *  it  *  *  k  kk  k  a  k  k  k k  kk  kk  kk 


3 12'U; 

3  I  30C 

3140  SUBROU  1'lMi  LI'iJ3 

3  14  20 

3  I  4  4C  NEGATIVE  HI  NON  I  A!.  APPROX  IMATION . 

3I46C 

314  1  COMMuN/ LPGPRM/ I  I  ,RLI  ,A1.U1  .  II)/,  108 

318!)  MMTF(o,m  "l  ilS  PROGRAM  COMPUTES  LPG  US I  JO  SCALED  BIN" 

3188  13  FORMAT! V) 

3101  WRITE! 0,13) 

31/0  WRITE!  6,  1  3  )  " TllFTA  =  ",  Tl  .  "LA  MliDA  ="  .  MU  ,  "ALPHA  =  '• ,  A  I  ,  "  BF  T  A  =",B1 
318  )  WRITE!  o,  U) 

3190  WW I TF ( 6 , 13)  "MFC  BINOMIAL  PROBABILITIES" 

3200  WRIlEt o,1 3) 

3210  WRITE! 0,13)  "  X  P  F" 

3220  WRITE!  6, 13) 

3230  C  =  1 / ! I . -T 1 ) 

3240  RKI  =  -T1*RLI/AL!X3!  I-TI  ) 

3250  C4  *  ( B1/! RKI +81  )  )**A 1 
3260  W  =  RK  I  /(  RK  I  +B  I ) 

3270  C2  =  0 

3280  DO  3850  MM=I  ,1001 

3282  M=MM- I 

3290  PI  =  1 

3300  IF!  M.NE.O  >GO  To  33/0 
3310  Pi  =  C4 
3320  A  =  - 1 

3330  B  =  I  FIX  (C/2  +0.5) 

3340  S?  =  PI/(C/2+1 ) 

3380  GOTO  3440 
3360C 

33o5  33/0  CONTINUE 
33/0  IX)  3390  1 1  =  I  ,  M 

3375  I=II-I 

3380  PI  *  P 1 ★ ( ( A 1 +M- I— I )  / !  M— I )*W) 

3390  3390  CONTINUE 
3400  PI  *  PI*C4 

3410  IA  =  IFIX!(2*M-I )*C/2+0.5) 

3411  PA3I A 

3420  IB  *  IFIX((2*M+I )*C/2+0.5) 

3422  Rb»IB 

3430  S2  •  PI/C 
3435  3440  CONTINUE 
3440  PL  3  *  C2 

3450  DO  3520  !X«IA+I,IB 

3455  X»IX 

3460  F3  *  C2  +  S2*(X-PA) 

3470  F2  «  F3  -  RL3 
3480  PL 3  »  F3 
349 OC 

3500  W.RITE<6, 1  3)  IX,F2,F3 
3510  WRITE  (8,13)  IX,F?,F3 
3520  3520  CONTINUE 
3530  C2  »  C2  ♦  PI 
3540  IF  C  C?  . OT .  .99  )0()  TO  3560 

3550  3850  CONTINUE 
3555  35  60  CONTINUE 

3560  WRITE  !8,13)  -99, -99,-99 ," END  OF  NB" 

35  70  RETURN 
3875  END 


LPGFOR.S 

LPG3 


353  )C**  **************************************** ***LP(] 4 

JWO  SUBROUTINE  LPG1 

3600C  LP:J4.S 

361 OC  COMPUTE  REOUIMTION  SI/I:  E  VS  THFTA  TABLF 

3620C 

3630  00  36/0  1=1,100 

36 35  T I = I  * , 0  I 

3040  H  =  -(  l-TI  )*ALOG<  l-TI  ) 

3tf> 0  R  *  11/ 13 

3600  WPITF<6, 13)  TI.R 

3665  13  FORMAT!  V) 

36/0  36/0  CONTINUE 
3680  RETURN 

3690C  ******************************************************* 

3/OOC 

3 /IOC  DEFINITIONS  IN  NAHM  I A  PROGRAM  LPGI.S 

3  12  OC 


3  730C 

LAMB 

RL  1 

3  7400 

ALPH 

A 1 

3  75  OC 

BET 

Bl 

3/OOC 

MEAN 

nun 

3  /  fOC 

CNST 

Cl 

3  /80C 

LIM 

RL2 

3/90C 

CUM 

C  2 

3800C 

FACT 

FI 

381  OC 

SUM 

SI 

3820C 

KN 

RKI 

3830C 

PROD 

PI 

3840C 

CN 

C3 

385 OC 

CONST 

C4 

3860C 

LAST 

L3 

3870C 

FX 

F2 

3880C 

FFX 

F3 

3890C 

SL 

S2 

3900  END 


LPGFOR.S 

LPG4 


CM  ( '  >/Y  I  i  \  JE  jC  rI  Rl*  1 1  u  = 


IHVr/ 


|  i;  J  =|  INVR/OiiJ/S  JHfOP.  ()(  ,1CM,  NO  JO) 

20-VH.R0P.S  —  CALCULATE  TCALED  NK.AlIVR  BINOMIAL  C  ‘Mil  CAL  JOINTS 

:v  x; 

4  )  SvJiiRODi  INI:  oNI.RoPlAVE-'FO.ni.XLT.GVI  r.HCW'l.I^JP) 

•ux: 

OO  COM;VON/UT/I  U(2o  > 

/MU 

otx:  >-;s rwM'r  lpc  parameters 

7>0C 

I0O  IWi  V)  =  I  W  i  (  >) 

I  IOC 

IPO  IKMro.Oi.  ))  UNITE  (o.l  3>AVEREO.()l  .ELT.CVLT.FCRIT 

I  JOO 

I  40 C  CHECK  FOR  NO  JLMATH  OR  NO  PENALTY  FOR  SHORTAGES  CASE 

I  hOC  1 N  THIN  EVENT.  >  <T  RE-ORDER  POINT  ft)  ZERO 

I  60  C 

I  70  I  F(l)1  .L>:.0.  >  OO  TO  601 

IK)  I  E  ( FCRIi'.LE . O.  )  GO  ft)  oOl 

I  90  C 

pax: 

Pin  I  E(  ELT.LE..  I  )RXT=.I 
220  IF<  f)l  .LE.O.  I  )  l»I -0.1 

230  IFICVLT.LE.o.ol  )CVLT=0.OI 

240  CALL  BSRCHf  AVHRHQ.Q) 

250  T  I  =  G 

2oO  S=CVLT*FLT 

270  B»S*S/ELT 

2hc)  ui=pi*r:lt 

290C 

300  Blal./B 

310  AI»»l*FLT 

320  RI=-TI/({ l.-Tl >*ALOG( l.-Tl >) 

330  RLI«U!/<ELT*RI  > 

340 C 

350  IF(  IWT6.GT.0)  WRITF(  6. 1 3)  "  Bl  Al*. 

300A  *  RI  :?L1  C  RKIM 

3/0  13  FORMAT!  V) 

3H0C 

390C  ESTIMATE  SNB  PARAMETERS 

400 C 

410  C»l. /C l.-Tl) 

420  RK  I  »-Tl  ★RL  I  /( ALOG(  l.-Tl)) 

430  IF(INT6.GT.'))  WRlTETo.o  J)BI  ,A  I  ,RI  f  RL  » .C. RKI 

440  63  FORMA  T (  8F9. 3  ) 

4  50C 

4  00C  <PCX*0. >-P0 

4  /OC 

400  P0»<B1/<RKI*BI  >  >**AI 

490  W*RKI  /<  RKl*Bl  ) 

500C 


SNBROP.  S 


S'  r  UP  IMTIAL  PROBABILITIES 


b  IOC 

s  c<  »c 

S  ♦)  x=o. 

b40  Xl.=0. 

->bO  FXL=0. 

*i  CX)  PX=PO 

b/O  IP(PX.LI V).  )  GO  TO  /oi 

i)Wi  xn=c/p. 

b/.)  FXN  =P() 

O(X)0 

oioc  is  cum >fcrit 

6POC  Ir  SO,  EXIT  L  >OP. 

OJOC 
O  400 
6bOC 

OCO  IF<  IWTo.OT.  9>WRI  I  Mo,  I  3  )  "  X  PXM, 

6/Oi  "  FX  XL  F  XL" 

OHOC 

690  2B0  CON  r  HUH 

700  I  F(  (  IvUTo.Oi  .0).ANI).  (X.LE.PO.  )  >W  ?ITF(  6, 9  3)  XN ,  FX,  FXN ,  XL  ,F  XL 

/ 10  93  FORMAT! ) X,2(F IO,2,2FIO.S) ) 

7;>0  I  F(  FXN.  OF.  PC  PIT)  GO  1*0  410 

/  30  I F  (FXN.GT.  .999)00  To  b  10 

7400 

750C  COMPUTE  NEXT  MEG.  BIN  PROS. 

7  00  0 

me 

/HOC  GEUEEUE  SNtJ  PROB  TILL  X* 

/ 90C  interval  IS  FOUND 

800  C 

8  tOC 

MFD  X  *  X+I. 

d  JO  XL=XN 

B40  FXL=FXN 

880  XN=XN+C 

BOO  PX  =  PX  *  <  ( A I  ♦  X-l.)/  X  )  *  W 

8/0  IF(  PX.LE.O.  )  00  TO  70  I 

ddO  FXN=FXN*PX 

090  GO  TO  ?.?.0 

9000 

9  IOC  INTERPOLATE  TO  FIND  X* 

9  POC 

9 JO  410  CONTINUE 

940  XCRIT»XL*(FCRIT-FXL >* < XN-XL )/(FXN  -FXL) 

9S0C 

9  60 C  ROUND  TO  GET  FINAL  VALUE 

v  iOC 

9U0  IRoPMXCRlT+.b) 

990  RETURN 

I  OOOC 

IOIOC  HIGH  PROTECTION  LIMIT.  SET  ROP  To  UPPER  LIMIT  OF  X. 

IOPOC 

1030  b 1 0  I ROP* ( XN+O. b ) 

1040  RETURN 

I060C 

I060C  DEMAND  Off  FCRIT  IS  ZERO.  SET  R0P*0 

I0/0C 

lOrtO  601  CONTINUE 


SNBROP. S 


IKK)  IR0P=0 

mo  RFTURIi 

1  I20C 

II30C  PS  =0,  bO  DFFAIlLi'  CALC UFA  H  ON 

I  HOC 

I  ISO  701  CONTI NUI: 

1160  I  ROP  -  AVFRF0*D1*FL  L*(  I  .-M.*CVLT> 

I  I/O  WI1  r»*  <  6,  />  3)1  OOP 

I  ID  /23  FORMAT!"  ********SNHPOP.  PX=0,  SSO  SHT  R0P=",I3) 

I  IVOC 

I  200  RETURN 

1210  FN  ) 

I  220C 

1230  SUBROUTINE"  USRCIKRO.Q) 

I  24  OC 

I  260  COMMON/ I WT/I WT(20) 

I2o0  Il)/*I  WT(O) 

12/0  IP8=1WT(6> 

I  26 OC  - A  VP  RFQ  SIZE  f-QUA  I  I  ON 

1200C - FUNCTION  FNR(O)  ■  ~Q/(  (  1-0)  *Al3X>  (  I  -0 )  ) 

I  300C 

I  31  OC  3  FF  >1*12-2  WO  FOR  FNR(O)  OFF  IN  IT  ION. 

I  32  OC 
I  330C 

I  34  OC  - 

I  35 OC  BINARY  SEA  RCH  ROUTINF. 

I 360C  - 

I  3/OC 

I380C  SOLVE  FOR  THE  VALUE  OF  0  WHICH  OIVFS  AND 

I  3s/ OC  AVR.  RF01I  ISITION  SIZFOF  RO 

I  400C 

1  41  OC  SFT  112  FN!)  POINTS  FOR  SEARCH 

I  ^20C 

1430  K=  0 

1440  CP=.  969 

1 450  R9*FNR<Q>) 

1460  Ql *.001 
|47()  RI*FNR(QI  ) 

I 48  )  GOTO  2080 


SNBROP . S 


1  490C 

15PX  CHECK  IK  RO  . ! . f; .  R 

I  5 1  OC 

1520  1940  CONTINUE 

1530  IK  Ro  .or.  R  )C()  l'o  2nno 

1  54  OC 

155)0  RE-SET  TOP  OF  INTFRVAL 

1560  R9  =  R 

15/)  09*0 

I  580  GOTO  2080 

I  5900 

1  600  2000  CONTINUE. 

I  61  OC 

I  62 OC  RR>bET  BOTTOM  OF  INTFRVAL 

1630  RI»R 

1640  Q 1  =Q 

1650  GOTO  2080 

1  660C 

16/OC  HALVR  THF  INTFRVAL,  AND  RR-RV ALUATF  THR  FUNCTION 

I  6800 

](f?0  2080  CONTI  MUR 
1/00  K=K-H 

I  /I  0  Q=(QI+09)/2 

1  720  R=FNR(  0) 

1  730C 

I  74 OC  I  F(  R  IS  WITHIN  .001  OF  RO,  >G0  TO  RRTURN 

1  7500 

1760  IFUD8.0T.0)  W R I TR (6,23)  K,R0,R,0 

1770  2  3  FORMAT <4'  ***BSRCH~  K,  RO  ,  R,  0",  15 , 3F8.3  ) 

1780  IF <  ABS(R-RO)  .LT.  .001  )G0  TO  2180 

1790  I FC  K  .GT.25  >00  TO  2180 

1800  GOTO  1940 

1810  2180  CONTINUE 

1820  RRTURN 

1 830  END 

1840  FUNCTION  FNR(O) 

1850  FNR— 0/<<  1-C)*ALOG(  1-0)  ) 

1 860  RETURN 

1 870  END 


SNBROP  .  S 


;ai  ai  u<;/f  n  »•  o'-bcM  p r  m  i=  nnxi/m-MPk  j 


I  o*  Ji,‘,i'4=  JllMiC/OiJ.I/DI-MPU  <.m  ilCU,  Mono  ) 

,?<)*!  *,*3.  b 

50  SUBROUTINE  OFMPARt  li)pFR,  1  P4, I  PS) 

C  >(J  ****  **  *******  *  **  ************  ******************  ******  **  ******  ******** 

bO  PA  RAM  FIT")  Nf'OO-MU 

OOC  ***  *  **  **  *  **  k  ****  k  *  ****  *  ******************  *  *  ***  ****  ****************** 

to  common/ i g.ju  ;/iguug 

pc  (’owuri/mr  /i  I  i  *.*,;• 

9 0  COMMON  / 1  Tw^'-K/l’i  W!  l  K 

10  >  COMMON/ M  MM  ril/l  TMM I  I! 

1 1.)  common /i  ro  ni/iTcTk 

»  2)  COMM  ON /Jl  :i  I-Ak/I  1  YIAR 

I  30  COMMON/NHOl  NY/Nt-NlMV 

I  40  COMMON/NF  MM  AX  /NF  F  MA  X 

1  ■■>()  COMMON/NIP-M/NI  i'FM 

I  00  COMMON /PRI  2F  /PRI  ?F 

I  /O  COMMON/NIK:NT/Nl)RNT(  I  ) 

160  COMMON /I Ml-  4Nl)/Il)FMNi)<  I  ,NUCG) 

1  y()  COMMON/ lt<!i  riJR/IWt:TUR<  I  ,NOQQ) 

200  COMMON/ 1  NR  0/1  W0<  I  ,N000> 

2  10  If-  <  ICJB  IO.N'-.  I  )  GO  TO  2b 

220  IV=l  1IME/I  TYEAR  ♦! 

230  I  CT-  IT  1  Ml  /I TOT!/  ♦  I 

2  40  ITW=(  I  TlMH-(  IY-1  >*l  TYEARI/ITWEEK  +1 

?bo  write* 0,10  mY,n/r.iTw 

260  103  FORMAT! “  DFMPA  U - «,5X,"YEAR  NO.  ",  13,5  X . «OTR  NO.*\I3,5X, 

2/0* 4  "wl-i-K  N0.",I3) 

200  3b  CONTINUE 
2900 

300C  INCREMENT  PER I 00  COUNTERS 

3I0C 

320  II)PEM=I  DPI:  R+  I 

330  UO  100  N= I , N ITEM 

340  IQ1  Y-*Il)I-MNO(  N,  IDPER) 

350  I  RET*  I  RET U  P<  N,  I  OPE  R> 

300  I  RO“IREQ(N.  IIIPER) 

370  IFdGMJC.F  >.  I  )W.RITF:<6,80I3>N,N0ENT(N>,  IOPER,  IQTY,  I  RET,  I RQ 

380  8013  FORMAT!  4X, '****I)EMPAR~N«',  15,'  NUENT-',  15.  '  IQPER»',I5. 

390*  '  I  OFMND-',  lb,-'  I  RFTRN*' ,  15 ,  '  I  HE O*',  15) 

400  IK(NOPNT(N).LE.O>  GO  TO  100 

4  IOC 

420C  GEN c RATE  SERVICEABLE  RETURNS,  IF  ANY 

4  JOC 

440  IF<  I  RET.  LF.  0)00  To  30 

4bo  n-iriME+i  tmnth 

4 oO  CALL  ENTER < IT,4,N,IRET,0) 

470  J)  CONTINUE 

4  80C 

490C  REFLECT  REOUISTION  FOR  ITEM  N 

500C 

501 C 

502C  SET  REOUISTION  COUNTER  KNTREO  TO  ZERO 

503C 

DEMPR3.S 

DEMPAR 


'  ?  K  i  I  -  ■  ~  o 

'tic  IM  Ml  Y.i  !  t)  I  no 

|a.»  i=i) 

‘  I  V*  ) 

‘>-<)0  o.;  i  A  I  U  A  Rl  Oil  1  >1 11  )N  SIZE  I  R  13/  MONi'P  Cm  >10 


*  t/i  »•  : 
no ) 
n  /« I 
■i  O  i 
*.  >’ ) 
ui  K) 

o  It'.:: 
OCOG 
O  30  C 

•  J  400 

O  ">0 
C  CA.) 

^o)  i.; 
002 
O  /:)•; 

o<>O0 


BO  ii-Ul\:h  !|-M<  .  0  ) 

v:  a l  l  N'-'OnN  I  I  :  !.•' ■ ' 

1  i  LOo  =  I  Ml  *1  ...J 
It  =  M ) 

it  <  ]  rivMH.i  c.  i  oi  v )  ;o  m  (*' 

LIMIT  LAST  OFOUl  ■  I  !  I  Oli  .<)  1 1:  A  i  J  ( ’  i  A  l_  if'JlA‘J|)  IN  ON' 
•-UlAl.S  1C1Y. 

iMi  / -  1  1)1 
oj  CONI  L  ill  I- 

K  N  i  Kt:  J=  KNfREO  f  I 


It-  (NF'-MAX-SO)  REOUl  SI  IONS  HAVE  Hill 

GENERATE!',  NU  f  ALL  R'MAlMNn  1)1-  MANN  ON  THIN  REQUIMTON. 


/  00  c 
i  IOC 
/9.) 

I  30  C 
/AO 
/NO 
/  n(  i 

(  to  < 

/  MO 

/voc 
MOOC 
b  Ok; 
mZOC 
O.iC) 
c340 
<}‘)0 

I I  on  c 
h  /o<: 
a  hoc 

8900 

9000 

v  IOC 

vzoc 
9  jo 

V  40 
9  NO 

)00  HV2  1 
9  70 X 
9:30& 

9900 
t  OOUC 
10IOC 


IM  KNTWiO  .  LT.  ( NEEMA  X  -  NO)  >  Of)  To  )JO 

i  r  =  iory-i  roi 
WM  TEC  6,  /  DITIMF.IR,  IOTV 
/  1  FuRMAT(4X,'****l)FMPAR~  I  HMF='.  lb,  '  F.F.L 
'  OVERFLOW  REQ=',I‘>,'  OF  '.IN) 

,7t/rf»:(  o.O)l  3  )N,NDEMT  <N  ) ,  IDPER,  IQTY.IRET,  IRQ 

DFIF'JMINF  PRIORITY  OF  THI  S  Rl  Oil  I  SI  TI ON 
ASSUMING  SOS  ARE  PRIORITY  I 


H)  CONTINUE 
IPR=I 

IF(RANI>EM(  ,2).LE.  PRI2F) 


IPR-2 


DETERMINE  ARRIVAL  TIME  FOR  THIS  REQUISITION 

LIMIT  ARRIVAL  TIME  TO  NO  LATER  THAN  .5  WEFIKS 
I3I-FoR,:  Till-  END  OF  THE  QUARTER 

LM  =  lTrfTR-0.t>*inv  EHK 

IT-ITI IE+RANDFM ( . ? )*FLOA T( LM ) 

IF  <  IGHUG.E ) .  <  )  WRITE!  6,  8023)  IT,N,  R,  I  RO, 19,  IPR 

FoRMA  ff  4X,  '****OEMP AR — *  ITs/,1 10,  *  N*',I4,'  R*',FS.3, 

'  INITAL  RFQ-SIZE*',  14.  '  FINAL  REOS !  4. 

'  PRIORITY*',  14) 

PUT  REQUISITION  OF  F • E.L . 


DEMPR3.S 

DEMPAR 


I  i i  I'Ai  t.  1-N  ir  (  1 1  ,  I  ,N  ,  t  W.  li'  yi 

I  n  h)  T  1 n  I  =  [  1  J 1  i  t 

I  '  >‘i  J  1  J  (  I  I'M  r.J.  I  .  1  •  J i  r  >  O')  i'll  l>() 

I  <  si  >o 

\<>o  >c********Ct{,:.\  i  *  n ?•:>;»'  i  At.*  '-vi  -ji 

I  o/oc 

lot'.o  i  of)  ';oi,’i  iNij| 

if 'Vo  n  =  i  r  i  mi  -  i  totr 

IIOO  L‘AI.1.  HN[VO(  ir,  ),(>> 

Min  M-njIM 

I  120  If  Jl> 

1130;:  :ion  *  uri'Ri-o  lb  ut  i  om  i-n  u',H>. 

I  I4tx:  RA  if  !f- P,  NI-OBNI  ^NKWATFL  MHGUIN  RKO  SIZFS. 

I  lb 00 


DEMPR3.S 

DEMPAR 


Il<v>  ;  i  I  !  ;'■  51 1  >1,,0(  W,  JOKY,  I  R) 

11/  K'  Hi-.  IVUU'IP  A  l/l-'XJI  >  f  J  ION  SI/K 

fl CONi  l  M'llriDINi;  (O  A  C.jv  WNOH.  (t|-  N. 

!  »  >•)«; 


1  >();) 

DIM  -  t| 

!>I  \  1  UCA 

l 

(  •<) 

.0 

p.|M3<  10, 

t3) , 

II) OC  IO,M 

) 

1  21.) 

DAM 

1  UCA  i  /  1  ,  < 

« 

10, 

3  1 

,  loo,  u  6 

, 1 000, 

1162/ 

1  ^-'O 

1  A  l  A 

i:riA>.j/itu 

1 

•  « 

1  2  <  >\ 

.  4  ■  > ,  .  H  > 

♦ 

M*  1 

.  o 

♦ 

1  24  K 

.2  / ,  .V 

* 

.<V> 

t  • 

/  3,  .  80,  . 

21  . 

.22 

,  .98. 

« 

26, 

I 

.0, 

1  2S0< 

.  ’  1,  .5  1 

* 

.  M 

•  • 

M.  .K2,  . 

20. 

.22 

.  .2/, 

• 

22, 

\ 

.  00 

1  Po  )K 

.no,  .  3*< 

• 

.O') 

•  • 

.89 

,  M  3, 

• 

26, 

1 

.00 

12/.  K 

,  • )  1  8 

f  • 

'  i. 

.42 

•  • 

20  ,  . 

98, 

.  /O 

,  .  /  6 

.  .  22  , 

.  2/,  1  .Of), 

1  280H. 

.012 

f  • 

1 

.  n 

•  • 

3  / .  • 

4H, 

•  SS 

,  .  00 

.  .  /»  , 

.8/ , 1 ,O0, 

1  2  vo\ 

.  >1 

t  • 

O  1, 

.  (2 

«  * 

80,. 

60, 

.  f>6 

..><6 

,  *  26  , 

1  .OO,  1  .(X)/ 

1  300 

DAI  A  IN)./ 10 

**■  ? 

t 

1  3!  OX 

1  « 

7. 

1. 

/* 

1, 

1  32  OX 

’  « 

2, 

3. 

4. 

s 

.  <>, 

/ 

2, 

10. 

1  33  >X 

■  « 

4  « 

6, 

M, 

10 

IS 

,70, 

7*3. 

10, 

1  M'tK 

■  * 

s. 

IO, 

*  , 

30 

»3l>. 

30, 

40, 

/O,  1  00, 

1  3S0X 

*  • 

8. 

lo. 

l‘». 

20 

->s 

•  '  >  • 

30 

HO  , 

ISO, 

300, 

1  Jo  OX 

*  t 

*’ , 

10, 

•  *3. 

20 

'  1 

>• 

30 

,80, 

200, 

600, 

I  J/OA  I  ,  20,4  0,  0<  >,  80, 1  <>0,2  00,  J(X),  600,700/ 

i  3«oo 
I  32  <x; 

I400C  PSiAtll  lSH  DEMAND  PATFCATFOORY 

I  4  I  OC 

I  420  IK)  10  1=1,0 

1430  I M=I 

1440  M*IIJCAT(I> 


1450  IF(IOTY.LR.M)  (50  10  30 


1  400 

1  4/OC 

10 

CONI'  1 NUF 

I480C 

1  49  OC 

OK  TH  DM  INF  2F  OU  IS  IT  ION  SIZF 

1  S(K) 

30 

CONTI NUF 

1510 

DO  40  1=1 ,10 

1520 

IF(CPROH<  I,  IN). OF. P)  OO  TO  60 

1530 

40 

CONTINUE 

1540 

1  55  OC 

00 

CONTINUE 

1560C 

I5/0C 

SET  RHOU IS  IT  ION  SIZF 

1580 

I FI I.GT.i  H50  TO  80 

1590 

IP«1 

I  600 

RETURN 

1610 

80 

CONTINUE 

1  620 

IFII.LT, 10)00  TO  100 

1  630 

IR*IRG< 10, IM) 

1  640 

RETURN 

1  650 

too 

CONTI NUF 

I  660 

DPROB*CPR')B<I,  IM)-CPI}OB< r-l ,IM> 

1  670 

RPROB*R-CPR(M<  I -I  ,IM) 

1680 

DRQ-IRQU,  IM)-IRO<  I-I  ,1  M> 

1  690 

IR=IRO<  I-I  ,IM>«MFIX<RPR(>B(H>R<yDPR084-.5) 

1  700 

RETURN 

1  710 

FND 

DEMPR3.S 

GETREQ 


COMPU18  ORAND  -1  HAN  X9 


>  >  h  k: 
2  m-k: 

8) 

2  jo  > 

?  u  oc 
2  32  OC 

>  3 HOC 
■>  140 
c  )C 
2  1O0C 
2  J/OC 
2  380 
2  390 
.'400 
24  I  00 
2-120C 
2  4  300 
244  0 
2460 
2  4600 
24  /'OC 
24800 
2  400 
2800 
20  10 
2620 
28  30 
2640 
2660 
2  660C 
26 /(K: 
268 OC 
2800C 
2  600 
26]  0 
2  0204. 
2  030* 
2  640* 
2660* 
2660* 
2670 
2  6/6 
2  680* 
2660 
2  /10& 
2  /2  ).* 
2  73  >C 
2  /4O0 
2  76  OC 
2  /oO 
277  OC 
2  /HOC 
2/9  OC 


6o  cnijriN.il- 

xo=i  vc  9 

18  l  ll'M-  18  ON  1.7  I  OitSERV AXIOM,  ASSUME  R-.04  39 

18  !IPI  .O  I'.  I  )  OO  lo  69 

US'-'  AVf-V-VJE  99  VAI  III-', 

8  1  99  =  .  04  J9 

P9  =  <  X9-  I  )/ !<9 
09 =P >+ I 

18  P9  -(),  JS i-  AVE  I  1 8M  VALIIFS 

IF<I»9. 1,8.0.  )  (JO  M)  4  1 
00  TO  I  70 

COMPUTE  VAR<  X>  ESTIMATH 

6  9  CONTINUE 

V=0. 

DO  oO  I -I, MPT 
V*V+C( I)*(X! 1 )  -X9 ) **  2 
00  CONTINUE 

V-V/! FLOAT! MPT— I )) 

S-SORT!  V) 

IF  VARIANCE  >  (X9-1),  CONTINUE  TO  100.  OTHERWISE, 

A  NEGATIVE  BINOMIAL  MODEL  DOESN'T  FIT. 

I F  ( I WT  (  3>.  OF.  2>  WRITE!  6,  73 >NPT, X9, V ,S, S/X9 , ! X!  I  ),I»I  ,NPT) 

7  3  FORMAT (T20,  'NPf-',  T30,  110/ 

T20,  'MEAN-'  ,T 30, F  I  3.2/ 

T20, 'VAR* '.T30, FI  3.2/ 

T20,'STI)  DEV-', T30.FI  3.2/ 

T20, 'COEF  OF  V- T30, FI  3. 2/ 

T 20, 'X!  I  )«',T30,20F6.2J 
IF(V.0T.(X9-| ))C0  TO  IDO 
IF <  IWT!  3).GE.  I  ) 

WRITE!  6,  63INPT,X9,V,S,S/X9,(X!  I),  I-I.NPT) 

63  FORMAT!  4!' ★***'  ),'MEGHIN.  VAR!X  )  <»  MFAN  , 

T20,' NPT,  MEAN,  VaR, S  ID  DFV.C.OF.V  ^ , I  4, 4F l 3. 2/ 
T20,'X! I)*',20Fo.2) 

CHECK  FOR  CONSTANT  REO  SIZE.  IF  NOT  ,  CO  To  51 
18!  { V . GT •  0.06)  .OR.  !NPT.LT.4>  )  00  TO  61 
USE  CONSTANT  REO  SIZE 


DEMPR3.S 

NEGBIN.S 


■  »•:  >2 

1  y.  >  hi  !  i-  -i  -f  I;  :  }.,  lN<»  OS  ) 

1  2  -O  I  h  l  |  .10  ./■;)  OO  1  t:  l  tin 

3  •>>  ’  IF  (Fi'lF  .1  I  .0.11  >i  (i  U) 

33.1  K 


<  i  -  >c 

roil, i ! 

hw’-i.h: 

j  410 

I  34 

[0FQ3  =  I4  >iF  )S*i 

3  42  ) 

IMNRi-oS.OT.  I'X)) 

MFC 

3440 

I  I  'HO  J(  '4  Rl  0  )  )  =1  ♦  1 

1  IS  1 

PRM)J(NRFoS  )^l- 

J4O0 

1  ob 

C(’N  IT  N  U I: 

i4  /  iij 

If  (F.d'.'O.  >'*■*  ).Hi 

i  r 

(JO  To  yon 


]o  >no 


.3-1  >  nf 

iso  >  i  ?o  conti  nut 

4b )  • )  300  Com  I  NUT 

JOT'  -c 

ib.ioc  FIJI  LAST  POINT  INTO  PWEOS(I>  TABLE 

3b -1  )C 

3bbO  NPFQs=r4RF'T  3  +1 

.1)00  IFINRHOS.  IT.  1)0)  NREOS^IOO 

3b/o  IRFOoIRRIQJ  )=I  ♦  ! 

ibbO  P  FEQS(  N  RF.Oj)  =1  .0 

3b90G 

ObnOC  IF  IWT<3).t;E.2,  WRITE  PPObAHILIfy  ARRAY 

361  OC 

3o?0  210  IF  (  I  WT(  3).L1.  >;t;o  TO  2S0 

3630  W  P I  Tl:  (6,21  3) 

304.)  21  i  FORMAT! //T'-D, 'REQUISITION  Sl/.F  C.D.F.'// 
ioSOA  '  I  R  P(X<sH)'/) 

3oO0  1)0  2  30  I*1,NRFoS 

36/00 

3680  WRI  115(6,22  1)1,  IRR)S(  I),PREQS(P 
3  690  22  3  FORMAT (14,1  6,F  11.4) 

3/00  230  CONTINUE 
.3  713  2S0  CONTINUE 

J/20  IF  ( I WT  <  3 )  .OE*2  >WRITE<  6,  323) 

3/30  323  F0RMATC//T20,' INITIAL  REQUISITION  COUNTS'// 

3/40A  T/,  'QTR',TIS,  'UNH  S ',T2/, 'REQ'//) 

3  /SOC 

J760C  USE  MONTH  CARLO  To  SET  FIRST  NDHIS 

3//0C  RFQ  COUNTS  CONSISTENT  WITH  THE  NEG  BIM  ESTIMATES 

3  /W>C 

3790  DO  380  I « I, NDHIS 

3800  IOTY-IDEMNDU  ) 

3bl 0  KNT=0 

3820  I TOT-O 

3830  IF(  IQTY.  LE.O  )G0  To  3  /0 

384 OC 

3HbOC  GENERATE  REQUISITIONS  TELL  TOTAL  UNJTS«>IGTY. 

3860C 


DEMPR3.S 

NCGB1N.S 


TvT- 


rjyTo  . 


I 


$  ,'rV- 


C*. 


<V 

I 


2d(  1) 

2  81  > 
,-82'> 

2  tan 

2  84 1 K ' 
2  H50C 
2  880C 
2B/0C 
2880C 
2890 
2  VO  > 
vvi) 

2  920 

2930C 

2  940C 

2V50C 

2960C 

2V70C 

29H0C 

2  99  DC 

J000C 

30  IOC 

302  OC 

j030C 

304  OC 

3050 

3000 

30/0 

30dO<v 

J090C 


Ill.hO  =  ! 

PPF  »);.>!  I  )  =  1.00 
I  RIO-!  »  i  =  X9  ♦  .5 
GO  TO  3  10 


COMPUTE  'I  KG  AT  I  VE  BINOMIAL  PARAMETERS 


100 


<JV  =  V/(  XV- I  > 

PV  =  09—  I 

RV= ( XV-  l  )/KV 

R2  =  (X9-I  > **2/! V-X9+ I  ) 


L  (:  i 


I  to 
1/3 


IkKOSCDMTH  REGUISITION  SIZF 

PRFOS!  I  >=PR<>8  ABILITY  REQUISITION  SIZF  <=IRFQS!I) 
NRFDS=NUMUFR  OF  POINTS  IN  REQUISITION  SIZE  C.D.F.  APq 

R9fP9.09*NFOATWE  BINOMIAL  PARAMETERS 

COMPUTE  NEGATIVE  BINOMIAL  PROBABILITIES 
USING  R9  AND  Q9 


INITIALIZE  VARIABLES 

CONTINUE 

IF! IWT! J/.CE .2 /WRITE (6, I  /3/R9,P9,09 
FORMAT! //'NEC.  BIN  PARAMETERS'/ 

'  RV=  '  ,F  8  .  3,  '  P9*'  ,F8.  3,' 


0P=',F8.3> 


& 

3100 

Y*0. 

$r\ 

3110 

P=l /( Q9**R9 ) 

3 

rl 

iv 

3120 

F=P 

3  1  30 

AVRa  1-P 

■v 

3U0 

P8=P9/Q9 

pp 

3150 

IP EQS( 1 >=J 

r- 

3160 

PREOS!  1  )*F 

31  tO 

NREOS* 1 

3  180C 

| 

3190C 

32QOC 

COMPUTE  NEGATIVE  BINOMIAL  PROBABILITIES 

-5 

321  OC 

MAX  RE  Oil  SIT I ON  SIZE  ALLOWS)  -  1000  UNITS 

3220 

DO  190  1=1 ♦ 1000 

i#* 

3230 

IPRNT=I 

L*. 

3240 

IF(  I. GT. 20 HPRNT-0 

>  % 

3250 

IK  MODI  1 , 1  0 ) .  EO .  0 )  I P  RNT*  1 

4 

3260 

IF!  ( I WT (  J  ) . GE •  3  )•  AND.  (IPRNT.EO.  1  > /WRITE! 6, 1 83  >  I .  Y.P ,F ,  AVE 

Jy3 

32/0  183 

FORMAT!'  I.Y.P.F.AVE-',  I  10,  F  10. 0, 3F 1 0. 4) 

*  • 

3280 

Y«l-I 

3290 

P»  ! !  Y+R9 //! Y+l  //*P8*P 

3300 

F=F>P 

3  310 

AVE-AVE  ♦  ! I-F ) 

i 

332 OC 

1 

3330C 

RECORD  APPROXIMATION  TO.F(X)  AT  LEAST  EVERY  .01 

v-; 


Ldvv.av.- 


DEMPR3.S 

NEGBIN.S 


i 

** 

a 

i 


<r. 

'*  i 


it'/)  <•>()  l'=,,,Avl!)‘-.H  .  .W 

ow  >  ;:ai  i.  nm.'Iw  l  <  »\  I  u  ) 
i  v.' >  rmr=i  ror+i  w 

i  voi*  KN  l'=\N  i  *  i 

jvio  IF ( I  fo f.  Li.  1  Oi'Y  )(;<>  I'd  ISO 

3^00 

WJOC  SIT  NtlMiJFU  OF  IJITIONS  KOUAt  lo  KNT. 

.W4  0C 

WSO  WO  !'?'•••')(  1  )  =  “<NT 
.wooo 

j//)C  *<W|  If-  I jFHU(>  MFSSAOF. 

JVOOC 

O'O  IF  (  I  VT(  JI.OH.DWWI  [fe'(o,  1/1)1  ,  IDFMMK  I ) ,  I  WFQ(  I  ) 

4000  ,i  /J  FOOMAi  (  H  10) 

4  o  I  OC 

40^0  i.m  coiriNUF 
4  030  WH  O  UN 

4040  FND 


DEMPR3.S 

NEGBIN.S 


'/  I  |  I'j:-  Ni  i MU  I  (  W,  I  It  > 

-'si  V'  A  *  A  *  *  *  ».  A  x  A  A  A  A  A  AnxAAAAAAAAAAA  A  A  AAAAAxAAAAAaAA  AAAAAA  AAAAAAAAAAAAAaAAAAI 

•i : ;v i ••  /as  A-lPTFR  'I  JO 

‘Ilf'  )•.'  A*  AAAAAA  A  A  A  A  A  A  AAAAAA  AAAAAAAA  AA  AAAAAA  AA  A  A  A  A  A  A  AA  A  A  AA  AAAAAAAA  A  A  A  AAA  A  A 

•;  i »  h; 

i,\)  oo.v.m  ; / i ., !  / 1 fi r ( «*o j 

‘I  I  J>  COM  i;. /\ !•/  UJO/I  l>i:MN!)(NOOO> 

4  M  )  Cv  \V.\ 

‘I  ! 1  >  i )  PI  V- Ob  lu'j  <(  10), C(  10), P(  10) 

4100  0C..1M  >N'N»-  Mill)/  NIM-GG,  IIP:(ji(  I0O)  ,Pia-<)S(  I  oo) 

41/0/  P 1  Vl-f;  R,  A  11(0,1)  P^I-IIDO  RANDOM  NlJ'AiiM', 

■»hOi:  l|  IHRVf  NH  1  HI-  CORRFbPONUI  NO  RHOU I  >111  ON  blZH  I'?. 

4  I  y'OC 

4/00  1)04/  )  [=1,  JlO-Ob 

4/10  II  (R.I.H.P I  )  HJO  JO  4  0 

4//0  4/  )  CUN  1'lN.Jft 
4/JO  4  j : )  CON  rifitJI- 
4/4  )C 

4/50C  bio  inr  m-oui  >i  rioN  si/h  i it 

4  zcoc 

4  ?  / 0  If*  (1.01,1  )'  >U  IO  40  ) 

4/dO  I R=  I  RH  Jb(  I  ) 

4/90  Rf-TURN 

•}  JOOC 

4  J I  0  4 1>0  CON  IT  NUH 

4  J/0  II-  CI.LT.  NRIOSHJO  IO  4rt(> 

4  J  JO  I  M*I  RHOS(NRRQS> 

4J4  )  RTTURN 

4  Jb  OC 

4  300  4 JO  CONTI NUH 
4  3/0  I M I  *  I  -  I 

4  3d0  IRI=I RHOb( l ) 

4390  IRIM1*IHH</{(  I M I  ) 

4400  IRDIF*  I  RI-IRIMI 

4410  IHCI  RDIF.Ol.  I  )(J()  TO  510 

4420  I  It*  I  HI 

4430  Hf-TURN 

4440C 

4450C  I  NTHRPOLATR  TO  DFTRRM  INR  REQUISITION 

4460C 

4470  b I  0  CONTI NUH 

44W0  DP  ROB*  PR  HQ  i(  I  )-PREQS(IMI  ) 

4490  RPROU  =  R-PR':OS<  I  M  1  ) 

4S00  IR*IRIMIaTFIX((  RPROB*I  Rl)  IF  /DPMOB  )*().b  ) 

4b  10  RETURN 

4b?0C 

4b2b  FIJI) 


DEMPR3.S 

NEGBNl 


4630C 

4*540C 
4560 
46/0C 
4  580 
4590 
4600 
4605 
4  610 
4  020 
4  630 
4  640 
4660 
4  600 


- RANDEM— U<0,  I  )  RANDOM  NUMBER  GENERATOR. 

RANDEM  —  U<0,  I  )  RANDOM  NUMBER  GENERATOR. 

FUNCTION  RANDPM(X) 

A  CALL  WITH  X  <  0.  INITIALIZES  THF  RANDOM  NUMBER  STREAM 
IE  ( X  )  10,20,20 
20  HN=R.IO*RAN  )EM 
RN  l*AMOL)<  RN  ,  BN  ) 

RANDKM *RN I /BN 
RETURN 

10  RHO*  /  ,0**1  1 
BN  =  I0.0**10 
RANDEM=-X 
GO  TO  20 
END 


DEMPR3 . S 
RANDEM 


•JS.- 


— *  k  ”  »  '  k  *  »  1 


-V  -  *._* 


CAl  Al.  t  *  ;/F  11  I-  1 'l-.iCPl  » 1  r  in\i 


I !  ■  *  >  / 1  i  M  1  1  M .  A . 


1  I',*  *  I ' 

2  i )«  i  ;j 
>• » 

*?  ■  )c  ~ 
so 

c>!  i:  — 
/.) 

?•(' 
v  O 

1 1  n 
i  le 
I  >1! 

i  *> 

i  a>  * 

ISO 
I  <JC ) 

I  H\ 

I  HO 

1  VO 

200 

210 

220 

2  10 
2  <10 

•2  SO 
2U> 

2  /O 
2  HO 
2  90 
300 
310 
J20 
.330 
340 
3  SO 
360 
3/0 
380 
390 
400 
410 
420 
430 
440 
4S0 
460 
4/0 
480 
490 
...  S00 

-::-5io 

520 
S  30 
540 


w=jilMiu/(ii.i/!  ;  itm  *.oi  hci>,  m/d 
114I.S 

suii|/o  ii  i;j*  itiiiiM 


PAPAMf  I  !•!/  M  *  Hi 


( :  1 1 A  P  A  C  r 1  1  f  A 1 . 0 , 1  ‘ ’  J , :  i  M  f  •  •  ( )  i  j  N , ; :  r  t ; :  I 

COM,1. ‘ON /A  1.0,1-  r:(4  ) ,  IIM  ,  NOIJN  (  2  ), 
COMMON/  ::  •! 

COMM'  )'J  / 1 0!  :.lt  5  / 1 !  >i'»ll(  • 

COMMON /I  I'OAV  /Ml  AY 
COMMON/ I  !  M  l  ‘  il/1  I'MrU  II 
COMMON /Ml  I £M/C)  i PM 
C  OMMON/  M I  )I:M/M  I Jl-M 
COMMON/NUUIS/NDm I S 
COMMON/ 1  Ni .11/ INI.' I 
COMMON/ IN  TYPE /IN  1  YPE 
COMMON/I  EdllO/l  f-.JHO 
COMMON  / 1  HOF  / 1  BOPUi  K  i),  !IH)P0P(  3) 
COMM  >N/I  ME  AND/l  I'l  MNl)(  I  ,N000) 
COMMON/FMII'  OS/FMFI  0S(  I  ) 

COMMON/ I  .W  u.i/invaci'U  ) 

COMMON /NON  >r  1/  lOplllM'd  > 

COMMON /f!  OF  MAC /I  MM  AC  (  I  ) 

COMMON/  JMnV.C/ir-N-lACt  l> 

OOMMO  J/NFP  JAC/NNI  OAOI  I  ) 

COMMON/ N0FVMNI)/NI)1; MNIK  1 ,N000> 
COMMON /N  PI  TiH?  /NP  !-TUP(  I  «  NOoQ) 

CO MMOil/NPI  l/NPI  ti(  I ,  MOOT)  > 
COMMON/Ni)»MT/IJ.>hNT<  1  ) 

COMMON/1  NVIjI) •'/!  NVIRJE  (  I  ) 

COMMON /Nd! )F I  /NlIOFT (  I  ) 
COMMON/NHOril/NHOTlH  I  ) 

COMMON/  41301  tl/NBO  I  U(  1  > 
C()MM()N/NL1()IR/NH0IF(  I  ) 
COMMON/.NHOTP/NU')]  PI  I  ) 

COMMON/W  QSI  /./PFGNIZ<  I  ) 

COMMON  /REO  lAO/Ff  OMAlX  I  ) 

COMMON /LTPROfl/l. TP WOU<  I  ) 

COMMON /I  TAOM/l  I  ADM(  I  ) 
COMM'lN/JCOST/UCOSl  (  I  ) 

COMMON /Al)>?/Al>P<  I  ) 

common / 1 Sul/  i su i. <  i ) 

COMMON/ I  PEO  /INI  0<  I.NOOO) 
COMMON/ 1  PE  TIP  /I  NFT(JN(  I ,  NOQQ  ) 
COMMON/ I  l?L/ 1  FU  I  ) 
coMMsN/in./rn.<  i  > 

COMMON/I  POL/ 1  RO!.(  I  ) 

COMMON/ 1  POTY/IPOTY<  1  ) 

COMMON /PM T  JR  /PMTHR<  I  ) 

COMMON/ RMEAN/PMF AN (  I  > 

common/ptp<-:no/ptpenih  i  > 

COMMON /RMA!1  /PMAIX  I  > 
COMMON/PE.7SUM/l?fcP5>JM<  I  ) 


,:i-*V'oUN,:.‘.<:rc:' 

;;(4  ) ,  UM,NUiJN(  2  ) ,  MO'l  Cl)(  4  ),  IOM,  IOP,  I  FIN  ,  I  PPPP 


IN1TM3.S 


550  COMMON /KN  i  /KN  J(l) 

5  00  COMMON )FUO<  3) 

5  /()  COMMON /ORO  )F/0 -T)OF<  i  ) 

Son  C0MMO  1/  I  Cl)f: '  >R/(C|>H)tf 

590  c 

o.X)  KNTW  INI  r 4 1 

OI0C 

o?.u;  :i»  T  —  oh  ~»;:ploi)S  OP  J  A  !  A  I 

o  JOC 

040  II)PI:p=N.)HM 

oSoO 

6  60C  PE-AD  DKMAND  I  >  A  J  A  POP  IT;  :M  N  FROM  LOGICAL  UNIT  LR 

0  /0  L  P«  I  Ni.U 

oa>  I  KN T—  I  KN  r+N  1  1 1-  •* 

090  DO  100  N  =  I  ,  N  I  TF  M 

/GO  10  CONTI HU R 
/  IOC 

720 C  URAL)  IT  FM  HAT  A  INPUT  PLUM  Mi  l  Ul 

/30C 

/ 40  i f ( intyvh. i- ).«.»)  oo  ro  20 

750C 

/ok:  PFAI)  HCf)  INPUT 

/  /OO' 

/HOC 

790  PFAIHLR,  8000,  l-Ul>  =  POO)  ALO.I  5N,UM  , BOOST! N)  f  NOUN,  MOTCI),  I()M.  ION, 

MOO  A  LIAPM(N),LTPP0D(N>.  !  PPL  ,  I  .’PPM 

H  IOC 

820  PI  PPPi/=P’LOA  I"(  IPPPP)/IOO. 

030  I  F (  I l;JtJ(».hO.  I  )rt<? I  TF<  6,8010)  I KNT,  Al  C , FSN ,UM ,  UCOS T( N  ) , NOUN .  MGTCD, 

8 40  A  I  OH,  i  OR  »l.r  AI)M(N  ) ,  LI  PPOIKN  > ,  IPPL.PIPPPP 

i}*;0  I  F(  I  RiiUO * P0«  I  )//  H  TM  0,801  5)  (1,1  =  1,10) 

H60C 

M/OC  PFAO  nl-MANf),  Rf-rriJRNb,  AND  'H-O-PRFOURNCY 

HBOC 

890  RFAIXLP,  800  >)< IDFMNOCN, J  )  ,J=I  ,  IUPFP) 

9 'JO  IF(  IFiPiO.FO.  I  IMUITHK  O.BOPO)  ( I  DHMNiX  N ,  J  >  ,  J  =  I  ,IT)PFR> 

9  10  PRAIXLP.OOOO)  ( IORTUP<N,  J),J=I  ,  IDPFP) 

920  I F <  IF. <00.1-0.  I  )URITM6,H030)<  I  PRTUPIN,  J),J«I  ,IDPfc*P> 

9  30  RFAlXI.'f.M'irv) )  ( IRFOCN,  J) ,  J  =1  . IDPRP  ) 

940  1F<  IFUJO.FO.  I  )WRITF(6,rt04O)(  IPRO(N,  J),J  =  I  , IDPRR) 

950  00  TO  50 

9  OOC 

9  /OC  RRAI)  BINARY  DATA 

9  800 

990  ?.<)  CONTI NOR 


INITM3.S 


»<  )  I\:Ji  ~  IKN  i  *■  I 

01  )  Rl -Ai;(!.7JAI  O.FSN,  f  R,  UM,  UCO 1 <  N  ) .  NOUN  .  Ml  , 1C! pill,  107, 

(V  )<  1.1AI>  UN), I  IRRO/HN)  .  1  i’H  ,*.•!  pp,  R 

i)i' )  IE  <  Iki  JG.>-  •).  I  >«  m',:(  o,  loin)  !  KNI  ,  Al V,Em4,UM.U0OST(  N)  ,  NoUN , 

(M'K  M  ( >  I  v  >,  10:1,  107, U  AM  U  1  )  ,  I.  i'PROI  <(  M  )  ,  1  »•*•!..  71  PPPP 

)•>:)  M-  V)(l.  7)  Al  C,  F.o  I  ,  I  7,1  U  MN'f) 

no  )  im  iPm/'i.i--’.  i )  i  r^(  o.Hopox  idi  mnix «i ,  j  >,  ./=  i ,  ippen) 

)/  )  REAiH  !.•>)  Af  >  ,  |  SIJ  ,  1  7 ,  I  RE  Hi  R 

OM  )  If-  C  I  EH  1C.  I  ).  I  KMl'l  <0,80  30  )<  17/  11  WIN,  J  ),  1=  I  ,  IDPFR  ) 

',70  RE  A!  Hi.  .') Al.  :.I  SN,  II/,  I  WFO 

1 0  )  im  ieoug.eo.  i  )f.:n  r  -  ( o*ho4c><  ire.hn,  i  >,  j=  »,  ii>per> 

11  x;  W IV  f  IE  INCUT  '  >A  I'A  To  FILE  07 

I  J'X 


]  tOC 

I4()C 

100  bOOO 
lo)  15010 
» 70S 

100  so  I  s 

1  VO  01)20 

30)  80-30 
21 0  8040 
22  OC 

2  JO  SO 

24 1  )C 
F>OC 
260C 
2  70C 

280  CALL  M KGB IN 
290C 

100  100  CONTI N UK 

31)  RETURN 
32  OC 

330C  INITIALIZE  ITEM  ARRAYS  FOR  A  NEW  SIM.  RUN 

34  OC 

350  ENTRY  IN  IT  42 

3oOC 

370  DO  210  N“  1  * NITFM 

JtiOC 

39 OC  INDICATE  ITEM  HAS  NDHIS  PERIODS  OF  DEMAND  HISTORY 

400C 

410  NDENT(N)*  NDHIS 

420C 

430C  ZER«)  DEMAND  HISTORY  RECORDS 

440  NRETACC  N )  *.') 

450  NDEMAC (N)=0 

460  NREQACIN  )»0 

4  70C 


Forma  r(  v) 

FORMA  I  (  /Is, IX.A2, 1X,A2, A4, A6,A3, IX,  A2,E I  I .2, IX* A6,A4, 
T5  3,  >AI  ,A4,A2,2X,2I  7,3I6,F5.2> 

Forma  r<  r;>i ,  loi  10) 

FORMA  f  (  *  IPE'MNO  DEMAND/QTR'  ,  (  12  I  ,  IO  I  I  0)  > 

Forma r (  'iritur  RETURN/0TR',n2i,  innon 
forma  re ^ i red  reo  /oiR'.rrei ,  ioi  io>) 

CONTI  RUE 

DETERMINE  PARAMETERS  FOR  NEGATIVE  BINOM I AL 
REQUISITION  SIZE  GENERATION 


INITM3.S 


•  «.•  -w-  Cv’  s,aV^S'  \V.  .nVvV.SV«.-.-  .*7 


1 

i  »,*  i 

1  OA  '  OM-ANP  II  Jio 

c* 

1  1  ><  V 

S' 

!  ‘>0  ) 

KK  -  :l  this 

1  ‘.1  ) 

1*0  1  I~l,i*tIIS 

£  v:'- 

1  ‘3,* '» 

Nl  H  m  |i :(  N ,  ts.r.  )  Tlit-AINIICN,  1  ) 

1 

1  >3  ) 

NR-  IMI’CIm,  KO-!  ,1-  I  I IRC  N,  I  ) 

s 

1  S-  ) 

NM-oiii,  Kk  )  ^  I  M-  )l  N,  I  ) 

■  J* 

$ 

1  .so 

KK=rv<-  1 

- 

IS<  )  1 

CONTINUE 

C-1 

1  S/  OC 

1  srtoc 

s»-r  j.jv  •  j  imry  r *m-—  in  i 

1  lwo 

/o  lfiVoj,;(;i >  =o 

1  oOO 

NORUt'TC  N  )  --  > 

1  • 


;-( 

si 

$ 


hi 

*1 


p. 

I  OMC 

1  o.'OG 

SET  INITIAL  iiACMWDl  R  C'HtNTI  Ri 

& 

1  03  )C 

1  04l) 

Nlit) I'll  (  N  )=  ) 

i 

1  6S0 

NFO I .)(  N  )^n 

1  O0D 

NbuIRCN  )-0 

1  6  to 

NbOTNC  N  )=0 

Vi 

>1 

1  08  0 

NiJOHTC  N)  =  ) 

1  09()C 

>♦« 

1  /OOC 

COMPUTE  DEMAND  RAT  lb 

9 

1  /IOC 

tm 

1  /20 

NN=N 

• 

1  /  30 

CALL  F  O  H5  /  M  NN  > 

,  * 

- ", 

1  /40C 

1  /bOC 

SHT  INITIAL  INVENTORY  ON- 

/OOC 

n  o 
/ao 

/vOC 
60  )C 
810 
H20C 
830C 
64  OC 
6*30 
660 
8/0 
180 
HOC) 
VOO 
010 
V20 
030 
940 
960 
960 
9/0 
080 
990 


HI.  1  =  (  L i AJ>M  (  N  >+l.TPROO ( W  ) )  *AI)R  ( N )  / 1  2 . 
IfcVACT(N)  =  WI.T  +0.6 


210  CONTINUE 

INITIALIZE  DROSS  ON -HAND  AND  ON-ORDER  STATISTICS 


120 


i\V 


7.4 


I  10  CONTINUE 

DO  130  N«I,NITFM 
I  F(  NDFNTC  N  )  t  LT.O  )  JO  To  130 
I F < I N  VAC T ( N  ) * LE . 0 )  GO  TO 
IBOHOHI I >  =1  BOH OH <  I )  +  l 
I  BOPOi I ( 2  )  ■  I  BOHOMC  2 ) +I  NV  ACTC  N  ) 

Il30H0.II  3)  -1 13  OH  OH  C  3>+IFI  X (  UCOST ( N>*FL<)ATI  INVACT(N  )  >3 
I  20  CONTINUE 

IF(  INVIHF(N).I  F.o)  GO  TM  130 
IBOFORC  I  )=  frioHORC  I  )♦! 

I  liOPORC  2 )  =  I  BOHO R(  2  )  ♦  I NV  DUF(  N ) 

I BOHORC  3)  -  I  BOHORC  J  >  ♦  I  FI  X(  UCOST  <N)*FLi)  AT  (INVDUECN  ))) 

130  CONTINUE 
140  CONTINUE 
RETURN 
2  000  200  WHITHC  6,8200) 

2010  8200  FORMAT!  MU  ,  ////20X,  'END  OF  FILE  RFADINF . '///* 

/A  Oil!  <  tt  I  f  * 


2020  CALL  OUT 

2030  CALL  ONTCST 


INITM3.S 


k>; 

►y 

2040 

CALL  PLOTIf 

& 

20b0 

STOH 

i 

2060 

FNO 

•  *• 
o, 

L": 

1  jL^k^au  .ttAjf a«l 

*»run=;hldg/obj/inord3.o(ucd/NOGo> 

*  I  NOW  D3  .  S 

SUBROUTINE  INORD(N/IOH/INIOR/IBO/IROL/IRQTy) 

C********** *» ft*********************************. **«*****•».* 

PARAMETER  NQJQ=38 

C  *****  *  *************  **************************************** 

c 

C  THIS  SUBROUTINE  DETERMINES  THE  RECEIPT  TIME  FOR  ASSETS 

C  ON-ORDER  AT  THE  BEGINNING  OF  THE  SIMULATION.  IT  IS 

C  CALLED  AT  THE  BEGINNING  OF  EACH  REPLICATION  AND  ITS 

C  OUTPUT  IS  STORED  IN  TWO  TABLES  (IDUF  &  IOWU)  TO  BE  USED 

C  FOR  EACH  SHORTAGE  FACTOR  WITHIN  THE  REPLICATION.  WHEN 

C  THE  ON-ORDER  QUANTITY  IS  ONE  AND  A  HALF  TIMES  THE 

C  EOQ  OR  LESS/  THE  TOTAL  ON-ORDER  ASSETS  ARE  ASSUMED 

C  TO  BE  THE  RESULT  OF  ONE  PROCUREMENT  ACTION  AND 

C  ALL  DUE  IN  AT  ONE  TIME.  WHEN  ON-ORDER  ASSEST 

C  EXCEED  ONE  AND  A  HALF  EOQ'S/  THIS  SUBROUTINE  THEN 

C  COMPUTES  A  DJE-IN  DATE  FOR  EACH  REMAINING  EOQ  SO 

C  LONG  AS  THE  COMPUTED  DUE-IN  DATE  IS  GREATER  THAN 

C  IFIVE  (16/533).  WHEN  THE  DUE-IN  DATE  IS  LESS 

C  THAN  IFIVE/  THE  REMAINING  ON  ORDER  IS  SET  AS  DUE-IN 

C  AT  IFIVE. 

C 

COMMON/ I DBUG/I DBUG 
C  OMMON/ IEBUG/I E  DUG 
COMMON/NDHIS/ND  HIS 
C  OMMON/ NDEMND/ N DEMN D ( 1 / NQ QQ ) 

C  OMMON/ N RE TUR/ NRETUR(1/NQQQ} 

COMMON/ITDAY/I TDAY 
C  OMMON/ I TMNTH/  I TMNTH 
C  OMMON/ L  TP  ROD /L  TP  ROD  (I  ) 

C  OMMON/L  TADM/L  TA9M( 1 ) 

COMMON/  IDUE/IDUECNQQQ) 

COMMON/ I0RQ/10RQ(NQQQ> 

C  OMMON/  J  CTR/JC  TR 
C 

C  WRITE  DEBUG  MESSAGE 

C 

I F(IEBUG.EQ.I)  WRITE (6/901 ) I  OH /  INIOR/IBO/IROL/IRQTY 
901  FORMAT(  ”  •  .  INOR  D  .  S  .  .  I0H«"/I6z"  INI  OR  =  "z  I  6/M  IB0*"/I6/"  IROL 

&I6z"  I RQTY*"/ 16) 

C 

C 

C  ZERO  ON-OROER  ARRAYS 

C 

DO  5  I«1/NQ0a 
I  DUE  (  I)  “0 
I  ORQ  (  I)  *0 
5  CONTINUE 

C - 

C  LOGIC  TO  TEST  LEAD-TIME  SENSITIVITY 

C 
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% 

V 


sV 


10PT  =  1 

IF  <  IOPT.  NE.2)  30  TO  7 
C 

C  SET  ALL  ON-ORDER  STOCK  DUE-IN  ONE  LEAD-TIME  IN 

C  THE  FUTURE. 

C  THIS  IS  A  VERY  PESSIMISTIC  ASSUMPTION. 

C 

JC  T R *0 

IF (INIOR.LE.O)  RETURN 
JC  T R *  1 

I D  UE ( 1 >*  (  LT  ADM ( N)+LTPR0D(N)>*ITMNTH 
1 0  RQ  (  1  )=  INIOR 
RETURN 
7  CONTINUE 


C - - END  OF  SENSITIVITY  LOGIC 

C 

C  LET 

C 

C  ISTKOBJ*ST OCKING  OBJECTIVE 

C  NET  A$STS*\|  ET  AVAILABLE  STOCK 

C  I  F I V  E  =  CLOC K  TIME  5  DAYS  INTO  SIMULATION 

C  C  HE  C  K IOR*M  A  X  SIZE  FOR  ONE  REQUISITION 

C  IHOLDIOR=J ORKING  VAR.  FOR  REMAINING  ON-ORDER  STOCK 

C  NREOQSOO-MA  X  REQUISITION  ON-ORDER  (INTEGER) 


I  OR  *  I NI 0  R 

I S  TKOBJ  *  I ROL  + 1 RQTY 
NETASSTS*!  OH*  I  OR  -  I B  0 
I  F I VE*5* I T  DAY 
CHECK  10 R*FLOAT  (  IRQTY)*1 .5 
I  HOLDIOR*IOR 
C 

C  IF  NREOQSOO  >  NDHIS  *  SET  EQUAL  TO  NDHIS  TO  AVOID  ARRAY  OV 

C 

K  QT R* 1 

C  COMMUTE  MAXIMUM  NO  OF  outstanding  ORDERS 

NR£OQSOO»IOR/IRQTY  ♦  1 
IF(NREOQSOO.GT.  ;;DHIS  INREOQSOO*  N  i*  H I  S 
N  DMD  S*0 

DO  50  J*1#NREOQSOO 
DO  10  K*KQTR#8 

N  DMDS*NDMDS+ND EMND(N#K)-NR£TUR(N#K) 

I AVAIL*NETASST  S+NDMDS 

I  F(  IEBUG.EQ.1  )  WR  ITE<  6/905)  J  ,K/ NDMDS,  I  AVAR 
903  FORMATC*  ORDER  N0*”#12#"  QTR*”#I2#”  NET  DEMAND*”#  16# 

ft  *'  I  AVAIL*”#  I  6) 

I FC IAVAIL.GT.I STK0BJ)G0T0  20 
10  CONTINUE 
C 

CALL  1 NT  I  ME (N#  R  L  T ) 

XMNTH*(RLT-24.) 

IDUE< J)*XMNTH*F  LOAT< ITMNTH) 

IF (IDUE(J).LT.I FIVE) I  DUE (J)*l FIVE 
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GOTO  30 


C 

C  COMPUTE  RANDOM  TIME  WITHIN  QTR  FOR  ORDER  RECEIPT. 

C  SET  IOUE(  J ) =DUE- IN  TIME  FOR  ORDER  NO.J. 

C 

20  T=RANDU(.2> 

KQTR=K+ 1 
0  RDQTR  =  K 

CALL  I  NT  I  ME  (N*  R  L  T  ) 

XMNTH  =  RLT  ♦  3.*(T-0RDQTR) 

X IDUE=XMNTH*FLOAT ( ITMNTH) 

I  DUE  C  J)  =  IF  IX  <  X  I  DUE  ) 

C 

C  LIMIT  DUE-N  TIME  TO  AT  LEAST  5  DAYS  IN  FUTURE. 

C 

I F Cl  DUE  (  J  )  .GT.  I  FI  VE  IGOTO  25 
IDUE(J)=IFIVE 
GOTO  30 
C 

C  SET  IORQCJ )*QUANTITY  FOR  ORDER  J. 

C 

25  I F(IOR.LE.CHECKIOR)GOTO  30 
IORQ ( J) = IRUTY 
GO  TO  35 
30  CONTINUE 

I ORQ  <  J )  =  IOR 
35  CONTINUE 

IFCJ.EQ.  NDHIS  >I0RQ(J)=IH0LDI0R 

I F< IEBUG.EQ.1 )  WR ITE ( 6#  904 ) T  *  XMN  TH  * IDUE (J )  # IORQ (J ) * IHOLDIOR 
904  FORMATC"  T*",F10.8#"  XMNTH*"*  F 1 0. 6*"  I  DU  E  (  J  )  ="  ,  I  7  , 

4"  I  OR  Q  <  J  )  »"  *  17#"  ON-OROER  BEFORE  PLACING  THIS  ORDER  *",I7) 

C 

C  UPDATE  COUNTER  VARIABLES 

C 

I F( IHOLDIOR- IORQ(J) >200,100,40 
40  IHOLDIOR*IHOLD IOR-IORQ(J) 

I STK08J=ISTK03 J+IRQTY 
I  OR* I HOL  D I  OR 
50  CONTINUE 
100  JCTRaj 
RETURN 
C 

C  PRINT  ERROR  MESSAGE 

C 

200  CONTINUE  * 

WRITE (6/21 3 ) 

213  FORMATC"***  *  *  I NORD ERROR  —  "  ) 

I DBU  G  *1 
RETURN 

END 
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*  I  NT  I  ME  .S 
C 

SUBROUTINE  INT  I  ME  (  N#  RL  T> 

C  THIS  RO J  TINE  COMPUTES  THE  ACTUAL  LEADTIME/  R  L  T  # 

C  FOR  A  GIVEN  REPLENISHMENT  ORDER  FOR  ITEM  N, 

C 

C  DEFINITIONS 

C  CP(I)  *  P  R  0  9 (  ACT  LT  TO  PRED  LT  <*  R AT IO(I)  > 

C  RATI 0(1) =  RATIO  OF  ACTUAL  LEADTIME  TO  PREDIC  ED  LEADTIME 

C  ASSOCIATED  WITH  CLASS  I 

C 

COMMON/IWT/  I W T (  2  0) 

COMMON/LTPROD/L  TPROO (1  ) 

COMMON/L  T  ADM  /LTADM(I) 

COMMON/ 1 RNDLT/ I  RNDL  T 
COMMON/I T  IME/IT  IME 
C 

DIMENSION  CP(1 0 ) #RAT 10(10) 

C 

C  THE  CP(I)  OATA  IS  FOR  A  GAMMA  PDF  WITH  MEAN  =  1 

C  AND  COEF.  OF  VAR.  =  .353.  »  / / CKCphK 

C  HENCE#  ALPHA  =  8#  BETA  *  1/8. 

C  THIS  IS  THE  MEDIAN  COEF  OF  VAR.  FOR  THE  62  ITEMS 

C  REPORTED  IN  APPENDEX  D  OF  HAYYA<1980>. 

C 

DATA  CP/. 000,. 3  7  7#. 2 5 5/. 4  90,. 700# 

&  .844# .927#.969#.987#1.00/ 

DATA  RAT  10/0.3, . 5/. 7#  .9,  1.10# 

&  1.3#  1  .5#  1  .7#  1.9#  2.70/ 

C 

I D  BUG* IW  T ( 1 ) 

C 

C  IF  I RNDL  T  *  1#  MONTE  CARLO  TO  DETERMINE  LEADTIME. 

C  OTHERWISE#  SET  LEADTIME  *  PREDICTED  VALUE. 

C 

IF (IRNDLT.EQ.1)  GO  TO  17 
C 

C  SET  LEADTIME  TO  PREDICTED  VALUE 

C 

RLT  *  LTPROD(N)  ♦  LTADM(N) 

RETURN 

C 

C  OBTAIN  A  U<0#1>  RANDOM  NUMBER  R 

C 

17  CONTINUE 
R»RANDU< .2) 

C 

C  FIND  PROB.  CLASS  THAT  INCLUDES  R 

C 

DO  20  I* 2 #  1 0 
I  F  (R.LE  .  CPU)  ) GO  TO  40 
20  CONTINUE 
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40  CONTINUE 


C 

C  INTERPOLATE  FOR  EXACT  RATIO  V 

C 

IL  *1-1 

RDELT*RAT IO(I)-RATIO(IL> 

CPOELT*CP<I)-CP<IL) 

C 

V=  RAT  1 0( I L )  ♦  < RDELT/CPDELT)*(R-CP(IL>) 

C 

LE  ADTM*  LTPROO(N)  ♦  LTAOM(N) 

RL  T*V*LE A  DTM 
C 

IF  CI0BU6. EQ.1) WRITE (6^1 1 3) I T I  ME , N, R , L E A DTM,V, RLT 
113  FORMAT!"  — t- I  NT  I ME - 1 T I ME ■  I  8, ”  N=",I4,"  R*",F4.2, 

"  PLANNED  LEAOT I  ME  13/"  V»%F4.2, 

"  RANDOM  LEAOTIME  ( MNT HS ) ■  F 5 . 2 ) 

IF  (RLT.LT..1)  RLT*.1 
C 

RETURN 

END 
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INTIME 


*  HFDi/OiiJ/EM  <'i'.(t(uCD/f.060) 

*  ^  N  T  R  L  3  .  S 

S  U lift  OUT  INE  ENT  EPLS<N#IQTY#I»RI#JT  IME) 
r  THIS  ROUTINE  U  A  C  k  ORDERS  REUUIS  it  lews  FOR  1TIM  N 

C  I  OT  Y  =UUAMT  I  TY  PLACED  ON  HACKGRDCR 

C  I  PR  1=1  HIGH-PRIORITY  REQUISITION 

C  I  PR  1=2  OTHERWISE 

C  JTIM£=  CLOCK  TIME  RFQ  WAS  RECEIVED 

COMMON / 1 1 bU  G /  I  DOUG 
C  OMMON/  N  f)MAX  /  N  b M  A  X 
C  OMMON  / NLOCt-K/  N  I.  0  C  d  k 
C  OM^ON/NBOIU/N  >30IU<  1  ) 

C  OMMON/NBOIk/N.iOIR  (  1  ) 

C  OMMON/  NBOTR/NOOTRd) 

C  OMMON/ NdOTU/N  OOFU(  1  ) 

C  OMMON/  NB  OPT  /NbOPT  <  1  ) 

COMMON/  l  3  A  C  P  T  /  I  fc*  A  C  P  T  (  1  ) 

COMMON /I C AN CL  /  I  C ANCL (1  ) 

COMMON/ I D  F  S  N  d  /  IDFSNO(I) 

C OMMON/ I LOCrfK/  ILOCBK(I) 

COMMON/  IPRIOR/  IPRIOR(I) 

C  OMMON/ IQTYb/I  Q  T  Y  B  (  1  ) 

COMMON/I  T  iMfc/I T  IME 
C  OMMON/ I TMbAC/  I  TMUAC  <1 ) 

C  RESERVE  A  STORAGE  LOCATION  FOR  THIS  INFORMATION 

I  PT  =  IL0CBK  <N!  0  CHK  ) 

f 

c  IF  BO  FILE  IS  FULL#  CANCEL  THIS  REQUISITION 

C  OTHERWISE#  GO  TO  5  AND  RE  CORD  f  H IS  0  6 . 

IF  (NLOCUK.GE.1  )  GO  TO  5 
C  THERE’S  NO  ROOM.  CANCEL  IT. 

IF (IDGUG.GE.1)  WRITE (6  #7 )  I T I  ME # I UT Y , I  PR  I 
7  F  ORMAT  (4  X# '***  *EfJTER8 — I  T  IM  E  *  *  #  I  is#  *  •# 

5  'BO  FILE  IS  FULL.  CANCEL  '# 

6  '  REQ  FOR  '#18#'  UNITS#  PRI*'#I5) 

C 

CALL  CUM(  ICANCL  #  I Q T  V  #N) 

RETURN 

C 

5  CONTINUE 

C  UPDATE  BACKORDER  COUNTERS 

C 

NBOT  R (N ) =NfiOT  R  <  N ) ♦  1 
NBOTU(N)»NBOTJ(N)+IQTY 
I  FdPRI.NE.I  )  o  0  TO  10 
N  80  I R (N ) =NbOIR  <N)+1 
N30lU(N)=Nb01J (N)+IOTY 
10  CONTINUE 

I F  < I OBUG.NE. 1 )  GO  TO  15 

WRITE (6# 13>N#I OTY# IPRI#NBOIU(N>  #N0OTU(N>  #NOOIR<  N>  #NS0TR(N>  # I PT 
13  FORMAT (4 X# *****  ENTERS  —  N«'#l5#'  IQTV«,#I5#'  IPRI=*#I5# 
ft'  NB0IU»'#I5#'  NB0TU»'#I5#'  NBOIR*'#IS#'  NB0TR='#I5# 
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1PT  =  •# I  5) 

IS  CONTINUE 

C  ■  DID  THIS  REQUISITION  CAUSf  THE  HACKORDFR  FILE  TO  OVERFLOW 

NLOCUK=NLOCti<-  1 
I F (NLOCbK.OE . 1 ) GO  TO  2G 
C  WRITE  ERROR  MESSAGE 

WRITE (6/91 ) 

91  FORMAT < 1HC#2D X#"ERROR — BACKORDER  fILE  OVERFLOW #F ILE  DUMP  ON"# 

8  "  NEXT  PAGF") 

WRITE(6#32) 

82  FORMAT( 1H1/1DX #  2  3  H*  *B  A  C  K 0  R  D  F  R  FILE  DUMP**) 

DO  83  K=1#NDMA  X 

83  WRITE (6 #84)  K#  I TMBA C < K ) # I D F S NB ( K ) # I  PR  I  OR < K ) # I Q T Yb < K ) # I B A C PT ( K ) 

84  FQRMATClH  #3X#7HREC  NO  =  # I  3 # 5 X # 7 H I T MB A C a # I  7# 5 X # 7 H I D F SNb= # II C # 5 X # 7 H I 
NPRIOR=#I1#5X#6HIQTVt»s#I7#5X#7HlBACPTs#I7) 

RETURN 

C  RECORD  QUANT  IT  Y#PRIORITY#FSN  I D#  AND  TIME  DATA  FOR  THIS  BO  REQ 

2C  I TMBAC( I PT)=JT IME 
I DFSNB( I PT)=N 
I  PRIORI IPTJalPRl 
1QTYB(IPT)=IQT  Y 

C  ARE  ANY  OTHER  BACKORDERS  OUTSTANDING  ON  ITEM  N 

I F(NBOPT (N).GT.O)  GO  TO  40 
C  RECORD  POINTER  DATA 

NBOPT  <N  )  =IPT 
IBACPTC IPT)=0 
RETURN 

C  IS  THE  NEW  30  A  PRIORITY  1  REQUISITION 

40  IF(IPRI.EQ.I)  GO  TO  60 

C  NOTE--  NEW  LOW  PRIORITY  BACKORDERS  ARE  INSERTED  LAST  ON  THE 

C  BACKORDER  CHAIN.  THE  REMAINING  STEPS  IN  THIS  PORTION 

C  OF  THE  SUBROUTINE  ACCOMPLISH  THIS  OBJECTIVE 

C  SET  JPT  EQUAL  TO  THE  FILE  LOCATION  NO  OF  THE  FIRST  BACKORDERED 

C  REQUISITION  IN  THE  CHAIN 

JPT=NBOPT(N) 

C  IS  JPT  THE  LAST  LINK  IN  THE  CHAIN 

49  I F ( IBAC PT < JPT) . EQ.O)  GO  TO  50 
KPT*IBACPT(JPT  ) 

JPT«KPT 
GO  TO  4  9 

C  RECORD  NEW  POINTERS 

50  IBACPTt  JPT)  =  IPT 
I BAC  PT< I PT )*3 
RETURN 

C  SET  JPT  EQUAL  TO  LOCATION  NO  OF  FIRST  BO  ON  CHAIN 

60  J  PT*NBOPT  <N) 

C  IS  8ACK0RDET  JPT  A  HIGH  PR  I  BO 

61  IF(IPRIORUPT)  .NE.1  )  GO  TO  80 

C  IS  JPT  THE  LAST  BACKORDER  ON  THE  CHAIN 

If (lBACPTt JPT) .EQ.O)  GO  TO  62 
K  PT* J  PT 

JPT*I8ACPT(KPT  ) 
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t,  0  T  J  61 

INSERT  NEW  3  0  AS  LAST  LINK  ON  CHAIN 
n2  I  0ACPT<  Jf'T)  =  F  T 
I UA  C  P  T(  I PT )=0 
RETURN 

IS  JPT  THf  ONLV  HO  ON  CHAIN 

80  I F  (  JPT.NE.N60P T <N> >  60  TO  81 

INSERT  NEW  30  AS  FIRST  LINK  ON  CHAIN 
NOOP T  (N  )  *  I  PT 
I HACPT(  I P  T )  =  J  J  f 
RETURN 

INSERT  NEW  JO  AS  LINK  BETWEEN  KPT  AND  JPT 

81  IBACPTt  K  P  T  )  =  I  3  T 
I 8ACPTC  I PT)=JP  T 
RETURN 

END 
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ft 


1  '*  *  ti  t-fint>:A);U/l>MAIi  UO(  UCD.ll!  IH  <ANI)U. O,  R 1 

2  )*»  it!-00/O0J/l>|- v.i*  iM.iWhl-l)  I/|lw>2l)A  1'A ••!);}“  ,  R 

• .  ■*'  'A  i  A  i>  <  .  S 

A  t  >. ' 

HUS  .VOtJlIN'  M-lFCIS  IN*.  >  l V-  WCORD:,  WII.J  AV»-RA<  MANUAL 

in  man. r;  in  on,-,  i-u  tha  r  s  if  in  imp  ran;f  rlu  in  run. 

! '  'C  IW'III'.]  I'll)',  r.UiNi:,  Air-  •  S  N--  RA  ri- ' )  IJV  MONO:  CARLO  to  UK 

mx:  :oNi»i;>n  n i'  ,m i  in  iinit  demand  rfooRUN. 


O'  >C 
O'  X  ’ 

/  •  k; 
•>tx: 

one 

i  .me 

I  IOC 

I  pne 
l  .y 
I  40 
l  oo 
i  t() 

I  to 
l  •/. 

1  VO 

2  00 
2  10 
?2f)C 
2  IOC 
240C 
200 
2  tOC 


nm;  NouiiNi-  riioir  pecorrd  a*4-  imnary  format  as 

DM- 1  OK.)  IN  I  N>>  I M  YO|  .  f,  Al’I'l-NDIX  A. 

2 A  NAM!-  Lp  0  NO  )Q  =.<  i 

1)1  MENS I OFsT  l-'.N  (  4  ),  NOUN(  2  ) ,  l  DEM*  U0(jU  ) ,  1 1?  1-1  (N  000)  .  I  RK0(  NOOO> ,  MOTCD(  4 
CHARACTER  MO  i  CD 
('AI  L  FPA  RAM  (  I  ,  I  J2> 

Rf:AI 1  (0,  i )  1.0  IT.LPRNl  .LSKIN 
J  FORMAT* V) 

KIM  NT,  "  LOUT  * ".LOUT ,  "  LPRNL  =",  LPf?MT ,  **  I.HKIP  =J,,LSKIP 
RFAIMS,  3)RLII,  PUB 

PI,’ I M  f ,  "  LOWER  HOUND  =  "  ,  RLU  ,  "  UPPER  HOUND  =",RUU 

INITIALIZE  RANDOM  NUMBER  SEED 
P— IMNl  JU(  —  .  I  ) 


*  ^ 

2  tOC 

IN  II  IALI/.F  COUNTERS 

X 

2  HOC 

5c- 

2  VO 

KNT«0 

300 

I0U1-0 

3  IOC 

n 

320 

ISKIP=0 

B? 

330 

10  CONTINUE 

E> 

340 

1SKIP  =  I  SKIP  ♦  I 

$ 

350 

IS  CONTINUE 

& 

360 

KNT  =  KMT  ♦  1 

& 

3  70 

IF  (MODIKNT,  10)  JtO.O) 

300  f 

3R0A 
400  I  I 
4  104 
420  23 

4  304 
440  JJ 
4  ,-»0  43 

460  S  J 
4  to  03* 
430  ^ 


rn  >  . 

L  *  .  -O 


RFAD(H,  END®  200  )AI  C,ESN,  IR ,U/4,  UCOST,  NOUN. MGTCD,  I  OH,  I  OR, 
LTADM,ll'PRO(>,  IPP.RIPP 

FORM  AT  (A  2,  A2,  A4,  A6,A3,  I  I  ,A2,F*>  .2,  A6,  A4 ,2AI  .A4.A2. 

21/  ,  3 1 2  ,  F  4‘.  2  ) 

I  FORMAT*!  /  ,  I  X,  ?A3,  A4,  I X  ,Ao.  A3,  I  X,  12,  A  3,  FIO.  2,  I X  .  A6,  A4  ,  I  X , 
2AI  ,A4,  A2.2I  /  ,2X,  313,  FS. 2  ) 

I  FORMAT ( A? ,  A2«  A4.AO.A3,  1 1  ,241  l) 
l  FORMAT* 3</IOX,UI  /),T9S,I  I0.TI2O.F  I  1 .0) 

I  FORMAT ( j(  /  (OX  ,  HI  /)  ,T  105,  110) 
k  FORMA  T(  3<  /IOX,  HI  /  )  .Tl  I  *> •  1 1 0  > 

hf:ai.(8  )ai.c,fsn,  ir,  idem 
RHAIHH  )AI.C  ,I:SM,  I  R,  IRFT 
RE  AIK  3  )ALC,FiiN,  IR.IREQ 
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••.vi  : 

*  4'u: 

•>•« ) 

•  >  a 1 

->  /'  i ; ; 

*-  >C 
va  ;0 
o  ji  >C 
o  I  •> 

<>.'0 
■>  <) 

6  40 

(>  Jfl 

060  6 
o  A  i 

ouoc 

6VOC 
/GO 
/  lo 
/  /OC 

/  *) 

/  40 
/SO 
/on 
/  A  i 
/t» ) 

ton 

hog 

H10 

H20C 

I  30  110 

640 J 

HbOC 

8O0C 

6  /OC 

800C 

BVOC 

900C 

9  IOC 

920 

9  30 

9  40 

950C 

900 

9  t() 

9  80 
990 

i  oooc 

IOIOC 
1 02  OC 
1030 
I  04  OC 


I  Fl 
I  Ft 


M 


•  >)■:  li 

U  I  * 


go i  *(  1 1 .  > 
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4  JlCIM  v)  .1-0.  “A, >.)/_«)  .()  l;t  I-, 


1)9 


4  >  )/> 


Sf-u-cr  i  ih,  .oi  im  ini  oa?j  ;r*  ii— i 


Il'*r 

[!.’=(. 

DO  HO  1=1,1 
ID=!IHI  IX-M(  1  * 

I  9  =  1  <>lfO-T(  I  ) 

1  COM  1  I  ill  - 

TI)  =  UC()VI>l-i  OA  I  (  1.1-19)/.*, 


IF<  ID. LI.  91. ,.»)»;()  lo  I*) 
II- ni).0i',  KOI-.  >  .,u  lo  1  '■) 


KI)FM=0 
K  PF'l  =0 
KRFf-O 

DO  no  1  =  I  ,M  )00 

<I)|:m=.<oFm>I  n'iff  ( i  > 

Kl?n  -KWT*  I  :>»•'!' <  1  ) 

KMFO^KNFU-M  IC-'Ol  I  ) 

I  K  (  Il)KM(  I) .  >T.  I  O. .»()())  do  TO  lb 
1F<  1  l?»~l  (  I  >.GT.  POOO)  GO  ro  lb 

COW  1'IMIF 

FBI  NT  THIS  itfm 


HFJ1-CT  TDll)  ITEM  IF  TOT Al.  DFMANDS  01)  RETURNS  EXCEED 
20  TIMES  THE  UI'FF R— rtdU Nl)  UfXJNJ)  KATE. 


RMAX=2>.  *IMh 

’?ATFaUGOST*FI.OAT  (  KOFM)  *4  ./FLOATS  NOOQ) 
IF ( RATF.GT. UMAX )  (X)  TO  15 

RATF=UCOST*FI.OAT (KRFT )  *4  ,  /FLOAT (  NOOO  > 
Ir< RATF.GT.  RMAX)  GO  TO  15 
IF (ISKIF.LT.  LSKIF  )  GO.  TO  10 
ISK 1 P=0 


GFNFRATF  REQUISITION  COUNTS  FOR  FIRST  8  QUA  RTF  RS 
DO  100  N= I  *  NOOO 
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L. 
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■ 
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I  ;  4  A.4 

)  (  1 01  Y )  IM  UP-  fJAWTE 

g 

1 1 1  >  x: 

1  HI  'b  ■  AftM.' 

r*i:  l,'l 

MMMTIOJ'.  bY  MoNi  -  CA 

is; 

i  i  ;>  >c 

$ 

M3) 

o'!  W=RAh>!.2> 

1  14  )(• 

i 

1  IS  ) 

CAi.i.  GFlR  -u(K.  i 

Ot  /  ,  1 

J) 

* 

II  O0 

i  kK  ii  -l  Kis •  J  i  H 

•'•J* 

I  1  V() 

nori.»i  i,)ii.  +i  (. 

1  lot 

li  ( i  ion.  .i  i . I03Y mo 

10  f>() 

‘.t; 

1  19;  k: 

*.♦ 

i  20'  k: 

i 

1  21) 

1  PRMH  »I  w<ill 

1  >20 

no  ro  1  »o 

E 

1  9  30 

Ml  I  RE  J!  1  J-O 

K 


&■ 


I  m 


>> 


94  0  I  JO  CON  Till  JE 

9t>0  Tr)Ol.=iJ(.\>Si*H  OAi'(KI)'-M) 

2ft O  W9ITM'>,2  JKfJT,  AI.C,l-bN,  l  R ,UM.  UCOST,  NOUN, MGi 'CD, 

>  /‘  )K  I  Oil ,  I  OR ,  LIAOM, LTPROO  ,  I PP ,  Rl  PIJ 


280 
290 
300 
31  J 
32 ) 

33  ) 

34  0 
JO  )C 
JOOC 
3/  >C 
33  0 
39  OC 
400C 
4  IOC 
420 


WRITKft.  /  i>  :  I >1 M.  K' l?f ■r.KRFQ, TOOL 

/  i  t'oi.*«nA  r<  uh,f  io.o) 

IK  IO'J  r  .01. I.PRNT  )(»()  10  200 
WRITE!  o,  //)  I  f)':'M 
WRITE!  o,  in  I  »?F-T 
IV PI  LET  <6,  /  / J  I  HFO 
//  FORMAT!  I  ox,  I  61  /) 


200 


CONTINWE 

OUTPUT  THIS  ITEM  TO  FILF  9 


v-. 

1  4302. 

i. TP  ROD,  I  l>P,  l< IFF 

1  440 

8  3  FORMA i'<  l  /,IX,A9,A2,A4 

% 

1  450* 

21/,  11  9,1-4.  2  ) 

51 

1  460 

93  FORMAT!  OX  ,81  /  ) 

1470 

WRI  TF ( V  JAM'  ,F"SN  ,9,1  DEM 

>  % 

1  480 

WRI  TE<  9  )ALC,F  »M  ,3,  I  RET 

:< 

1  490 

M RI TF <  9 )  A LC  ,E3N  ,4,1  REO 

$ 

1  500 

IOUT-IOiIT-H 

SI) 

S20 


S30 

*>40 

S50 

560 

0/0 

*>80 


300  CON1JN8F 

IF< IO JT.GT.LOUT>  c;o  TO  900 
GO  TO  10 
900  COOriHtJF 

write!  o,  903  jknt,  i  our 

903  FORMAT!  /////T40,' RECORDS  IM  «'  f  I5/T40,'  RECORDS  OUT  «'  ,  1 5> 
STOP 

f:nij 


DATAB3.S 


APPENDIX  B 


FORTRAN  SOURCE  LISTINGS 


o  o 


i. j > 


- 5FOIN  L ADR AND I  AN  LOOP 

DO  ?O0  MRUN»I  ,NRUN 
[RUN-MP1JN 

';OSHRT*CSHOR7(  MRUN ) 

I F ( IOBUC . FQ. I ) WRITF( 6, 8090  JCOSDRT 

N090  FORMAT!  //I0<  '*+**'),  'BFOIN  S  I  V.ULA  I  1 0*4  '  , 
V  VI  fH  C')‘>HRT  «',FI0.4.10< '★***')//> 

RFWIW)  ITFM  INRUT  FILF  INLN 
PFrtlND  INUJ 


- RFCOR0  RANDOM  NUMrtFR  SFFD 


—  INITIALIZF  RANDOM  NUMBFP  STRFAM 

INITIALIZF  DFMPP3  RANDOM  N.JMNFP  8  f  Ri  •  A  M 
IF<  f  I  SFFD.  f- 0.0)  RNLAST  -  tiANim\(-.?> 

IF(  nSFPD.NF.O)  RNLAST  *  !?ANDFM(-  ABS(  FI.O' T<|  [  ,  ) 

INITIALIZF  DFNFRAL  RAN'Vjm  NH.MHFP  SDMJ 


A/eu,'  hgdgv 
M4/AJ 


cepe 


RNLAST  =  1 1  SFFD 

IF  (  FI  SrFI)  ,f-0. 0)  RNLAST  =  RANDtJ(-.l) 

IF(  flSFFD.NH.O)  RNLAST  -  PANDIK -AUS(  RNLAS i  ) ) 


CALL  ZI;RO 


INVRSM.S 


